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ABSTRACT 

This paper presents the calibration of BVRIH-o s Tully-Fisher relations based on Cepheid distances to 2 1 galax- 
ies within 25 Mpc, and 23 clusters within 10,000 km s"'. These relations have been applied to several distant 
cluster surveys in order to derive a value for the Hubble constant, Hq, mainly concentrating on an I-band all-sky 
survey by Giovanelli and collaborators which consisted of total I magnitudes and 50% linewidth data for ^550 
galaxies in 16 clusters. For comparison, we also derive the values of Ho using surveys in B-band and V-band by 
Bothun and collaborators, and in H-bandby Aaronson and collaborators. Careful comparisons with various other 
databases from literature suggest that the H-band data, whose magnitudes are isophotal magnitudes extrapolated 
from aperture magnitudes rather than total magnitudes, are subject to systematic uncertainties. Taking a weighted 
average of the estimates of Hubble constants from four surveys, we obtain Ho = 71 ± 4 (random) ± 7 (systematic) 
We have also investigated how various systematic uncertainties affect the value of Ho such as the internal extinc- 
tion correction method used, Tully-Fisher slopes and shapes, a possible metallicity dependence of the Cepheid 
period-luminosity relation and cluster population incompleteness bias. 
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1. INTRODUCTION 

The Tully-Fisher relation, a luminosity-linewidth correla- 
tion for spiral galaxies, has become the most widely used 
method for determining galaxy distances independently of red- 
shifts. Although the concept had been used variously, for ex- 
ample, by Opik (1922) and Roberts (1969, 1975) over several 
decades, it was the groundbreaking paper by Tully & Fisher 
(1977: hereafter TF77) which laid out its potential as a pow- 
erful extra-galactic distance indicator The Tully-Fisher (TF) 
relation is an empirical correlation between the absolute mag- 
nitude of a late-type spiral galaxy and its maximum rotational 
velocity. The latter quantity is usually determined from radio 
observations of the neutral hydrogen 21cm emission spectrum 
or from the optical rotation curve and is distance-independent. 
The TF relation has been used not only to investigate the density 
and velocity fields of the large-scale distributions of galaxies, 
but also to measure the value of the Hubble constant. Readers 
are refeiTed to Dekel (1994), Strauss and WilHck (1995) and 
Giovanelli (1997) for extensive reviews on these subjects. 

After the first application of the TF relation using B-band 
photometry of spiral galaxies (TF77), it was soon discovered 
that there were mainly two advantages to using near-infrared 



over visible-wavelength photometry: (1) the uncertainties due 
to interstellar absorption in both our own Galaxy and in the ob- 
served galaxy would be greatly reduced; and (2) the infrared 
emission is dominated by late-type giants and thus there is no 
significant dependence of mass-to-light ratio on galaxy mor- 
phology. The first study to apply the H-band (1.65/im) - HI 
linewidth TF relation was reported by Aaronson, Huchra & 
Mould (1980), who used a large-aperture H-band magnitudes 
to obtain distances to the M81, Sculptor and MlOl groups. 
Aaronson and Mould (1983) later reported that in the infrared, 
the TF relation is morphology-independent, unlike in the blue 
where there is a small dependence. They also found that the 
slope in the relation is wavelength-dependent with its value in- 
creasing from the blue to IR. The scatter has been reported to 
be smaller for the infrared relation, making it a better distance 
indicator. Furthermore, the TF relation has been reported to be 
different for field and cluster galaxies or for galaxies in different 
clusters (Sandage and Tammann 1984). 

Major improvement in the TF relation has been achieved dur- 
ing the late 1980s using CCD images. As a result, numerous 
all-sky surveys of spiral galaxies were carried out which used 
the TF relation to estimate their distances, to probe the density 
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and velocity fields of the large-scale structure of the Universe. 
These include an r-band survey by Willick (1991) who used 
the 21cm linewidths to measure the rotation velocity, and also 
by Courteau (1992) who used the optical rotation curves. Han 
(1991) prsented an I-band TF all-sky survey, which concen- 
trated on the galaxies previously observed by Aaronson et al. 
(1982: hereafter AHM82) in their H-band survey. Most re- 
cently, Giovaneli et al. (1997: hereafter G97) presented an all- 
sky I-band survey of cluster galaxies out to 10,000 km s"'. 
Mathewson et al. (1995) also measured I-band TF distances 
to clusters in the Southern hemisphere. Although the data have 
not been published as of now, Willick (1999) has presented a 
peculiar-velocity study based on his Las Campanas/Palomar 
Survey which is comprised of R-band photometry and Ha ro- 
tation curve observationas of spiral galaxies out to 10,000 km 
s"'. Dale et al. (1999) also has presented the velocity field 
study based on their I-band spiral-galaxy survey reaching as 
far out as ^25,000 km s . Courteau (1999) has also recently 
completed V and I-band photometric and Ha rotation-curve 
observations of spiral galaxies. 

Although there is no firm physical basis for the TF relation, 
several groups have tried to recover the TF physical parame- 
ters using cosmological numerical simulations. Some exam- 
ples include that by Steinmetz and Navarro (1999) who used 
high-resolution standard CDM simulations. They reported that 
the slope of the TF relation can be reproduced in a hierarchical 
clustering model naturally, suggesting that the TF relation is 
a direct consequence of the cosmological equivalence between 
the mass and circular velocity, which was also the conclusion 
of Mo, Mao & White (1998). On the other hand, Eisenstein & 
Loeb (1996) argued that because the scatter due to variations 
in formation histories is larger than that observed, the TF rela- 
tion is a consequence of subsequent feedback processes which 
regulate the gas dynamics and star formation in these galaxies. 

The main objective of the Hubble Space Telescope Key 
Project on the Extragalactic Distance Scale (Hq Key Project) 
is to measure the value of the Hubble constant to an accuracy 
of 10%. This is achieved by first obtaining precise Cepheid 
distances to nearby spiral galaxies and then by calibrating sev- 
eral secondary distance indicators. The Hq Key Project has now 
completed the observations of all of its target galaxies and the 
Cepheid data and distances for each galaxy have been presented 
as part of this series in papers I through XXIII. This paper 
presents research and analysis which are part of our final effort: 
the determination of the Hubble constant through the calibra- 
tion of a variety of secondary distance indicators. In particular, 
this paper focuses on the calibration of the optical and near-IR 
TF relations using the Cepheid distances and the applications of 
the calibrated relations to distant cluster samples to determine 
the value of the Hubble constant. Other secondary distance in- 
dicators that are calibrated as part of the Hq Key Project are the 
surface brightness fluctuation method (Ferrarese et al. 2000a), 
Type la supernovae (Gibson et al. 2000) and fundamental plane 
(Kelson et al. 2000). The measurements of the Hubble constant 
will then be combined together in a consistent fashion in Mould 
et al. (2000). The extrapolation from a local value of the Hub- 
ble constant to the global value will be discussed in Freedman 
et al. (2000). 

We have obtained new BVRJ surface photometry of most of 
the galaxies suitable for the absolute calibration of the TF rela- 
tion. These data are presented in Macri et al.(2000) and sum- 
marized in Section 2. The photometric data will be combined 
with the rotational velocity data compiled from the literature to 



calibrate TF relations in these four optical wavelengths and in 
H band (Section]^). We also address the values of the intrinsic 
dispersion of the TF relations and the role played by a second 
parameter and non-linearities. The calibrated TF relations are 
then applied to distant clusters, mainly those from the I-band 
survey presented by Giovanelli et al. (1997a). The I-band TF 
distances to those clusters and values of Hq are presented in 
Section 4. They are then compared with values obtained by 
examining surveys at other wavelengths. In Section 5, vari- 
ous systematics that could affect the measurement of the Hq 
are examined. These include: (1) consistency between the cal- 
ibrators and cluster samples; (2) dependence on the TF slopes 
and shapes; (3) dependence on the internal extinction correc- 
tion method used; (4) metallicity dependence of the Cepheid 
period-luminosity relations; (5) large-scale velocity field of the 
Universe; and (6) cluster population incompleteness bias. Fi- 
nally in Section 6, our Hq values are compared with previously 
derived TF Hq estimates. 

2. DATA 
2.1. Cepheid Distances 

In Table 1, we list the spiral galaxies used in the TF analysis 
and the sources of their Cepheid distances. The list is com- 
prised of : (1) 15 galaxies which were observed as part of the 
Ho Key Project, (2) 3 galaxies whose Cepheid distances were 
measured by other groups using the HST, and (3) 3 galaxies 
with Cepheid distances determined from ground-based obser- 
vations. Note that there are several galaxies, such as NGC 4321 
(Ml 00) and NGC 5457 (MlOl), that have been observed by the 
Key Project but are not suitable for calibrating the TF relation 
because their inclinations are too small. 

The distances given in Table 1 were all taken directly from 
the references as indicated. However, the corresponding errors 
were re-estimated so that they are consistent from galaxy to 
galaxy, as the understanding of uncertainties in Key Project 
Cepheid distances have grown over the course of the project. 
The details of which and how uncertainties were treated are 
found in Ferrarese et al. (2000b). 

2.2. Photometry of Tully-Fisher Calibrators 

We summarize briefly in this section the surface photome- 
try observations and results, as this subject is presented in full 
detail by Macri et al. (2000). BVRI observations of the TuUy- 
Fisher calibrator galaxies were conducted between March 1994 
and May 1999 using the Fred L. Whipple Observatory 1.2-m 
telescope and the Mount Stromlo and Siding Springs Observa- 
tories 1-m telescope. Exponential disks were fitted to the sur- 
face brightness profiles, following Han (1992), to extrapolate 
total magnitudes. Both the profile fitting and the determination 
of the inclination were done interactively, in a manner similar 
to Giovanefli et al. (1997a). 

For the IR data, we use //-0 5 magnitudes presented by 
Aaronson et al. (1982), since with the exception of some work 
by Bernstein et al. (1994), that is all that is currently available 
in the literature, in terms of a consistent set of magnitudes be- 
tween the calibrators and cluster galaxies. 

2.2.1. Corrections Applied to the Photometric Data 

For the TF application, the photometric data were then cor- 
rected for Galactic extinction {Ag.\), internal extinction (A,„,.a) 
and k-correction. The Galactic extinction values were esti- 
mated using the lOO/^m maps reprocessed from IRAS/ISSA and 
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Table 1 
Cepheid Distances 





Tvne 


tn — (ni3§) 


Reference 


IT lUJv^ I 


NGC 224 


3 


24.44 ±0.10 


Freedman 1990 


ground-based 


NGC 598 


6 


24.64 ± 0.09 


Freedman 1990 


ground-based 


NGC 925 


7 


29.84 ± 0.08 


Silbermann et al. 1996 


Ho KP 


NGC 1365 


3 


31.39 ±0.20 


Silbermann et al. 1999 / Ferrarese et al. 2000 


Ho KP 


NGC 1425 


3 


31.81 ± 0.06 


Mould et al. 1999 / Ferrarese et al. 2000 


Ho KP 


NGC 2U9U 


5 


30.45 ± 0.08 


rnelps et al. 1998 


Ho KP 


NGC 2403 


6 


27.51 ± 0.24 


Madore & Freedman 1991 


ground-based 


NGC 2541 


6 


30.47 ± 0.08 


Ferrarese et al. 1998 


HoKP 


NGC 3031 


2 


27.80 ± 0.08 


Freedman et al. 1994 


HoKP 


NGC 3198 


5 


30.80 ± 0.06 


Kelson etal. 1999 


HoKP 


NGC 3319 


3 


30.78 ± 0.12 


Sakai etal. 1999 


HoKP 


NGC 3351 


3 


30.01 ± 0.08 


Graham etal. 1997 


HoKP 


NGC 3368 


2 


30.20 ± 0.10 


Tanvir et al. 1995/Gibson et al. 2000 


HST 


NGC 3621 


7 


29.13 ±0.11 


Rawson et al. 1997 


HoKP 


NGC 3627 


3 


30.06 ±0.17 


Saha et al. 1999/Gibson et al. 2000 


HST 


NGC 4414 


5 


31.41 ± O.IO 


Turner etal. 1998 


HoKP 


NGC 4535 


5 


31.10 ±0.07 


Macri et al. 1999 / Ferrarese et al. 2000 


HoKP 


NGC 4536 


4 


30.95 ± 0.08 


Saha et al. 1996/Gibson et al. 2000 


HST 


NGC 4548 


3 


31.04 ±0.23 


Graham et al. 1999 


HoKP 


NGC 4725 


2 


30.57 ± 0.08 


Gibson et al. 2000 / Ferrarese et al. 2000 


HoKP 


NGC 7331 


3 


30.89 ± 0.10 


Hughes et al. 1998 


HoKP 



COBE/DIRBE data by Schlegel et al. (1998). The /-band cos- 
mological k-term was adopted from Han (1992): ki = 0.l6z. 

The amount of internal extinction in spiral galaxies and its 
correction methods still remains a very controversial issue. 
Several papers on this issue include Holmberg (1958), Valentijn 
(1990), Han (1992), GiovanelU et al. (1994), Burstein, WilUck 
& Courteau (1994), TuUy et al. (1998). For the internal ex- 
tinction corrections, we have decided to apply those derived 
by TuUy et al. (1998: hereafter T98), which are functions of 
both major-to-minor axes ratio (a/b) and Unewidth, and are 
expressed as: 



, r98 



:7Alog(a/^), 



where 7a is a function of linewidth: 



7 = 



(1.57 + 2.75(log(iy20%)-2.5)) B 

(1.15 + 1.88(log(W2o%)-2.5)) R 

(0.92+ 1.63(log(W'20%)- 2.5)) I 

(0.22 + 0.40(log(W2o%)- 2.5)) K 



(1) 



(2) 



The idea behind the dependence on Unewidth is that the light 
needs to travel further on average in a larger galaxy, thus re- 
quiring a larger extinction correction. However, it is noted that 
careful examination by Wilhck et al. (1996) revealed no depen- 
dence of the internal extinction on linewidth or lumiosity. 

In Appendix A, we review two other expressions for internal 
extinction derived by Han (1992) and Giovanelli et al. (1994). 
The main difference in these two studies is that the extinction 
corrections are dependent on the morphological types of galax- 
ies. Mainly because the uncertainty in the morphological clas- 
sification of distant galaxies is large, we have decided to use the 
T98 formulae in this paper. Later in Section 5.7, we will exam- 
ine whether the value of Hq depends on the internal correction 
method. 

The T98 corrections in Equation 1 include no independent 
formulae for the V and H bands. For these, we adopt the ratio 
Ay = 1.5A/ and Ah = 0.5A/, which are reasonable values ob- 
tained by interpolating the T98 "extinction law" (Equations 1 
and 2). Note that these ratios are sUghtly smaller than those 
derived by Cardelh, Clayton and Mathis (1989). This is to be 



expected since Cardelli et al. used a fixed value for the opti- 
cal path length, when in reaUty, blue Ught is attenuated more 
quickly than red, and thus different optical path lengths are 
sampled at different wavelengths (shorter in the blue than in 
the red) (Han 1992, T98). 

Putting all the corrections together, the corrected apparent 
magnitude is then expressed as: 



ml = moh'i,\-A^ 



c.\- 



(3) 



Note that the essential point for present purposes is to adopt 
consistent corrections between the calibrator and the distance 
galaxy samples. We emphasize that the corrections do not need 
to be optimum; they need to be consistent. 

2.2.2. Inclination 

Accurate measurement of inchnation angles is crucial in the 
application of the TF relation since it is possibly the largest 
source of error in both magnitudes and linewidths. We mea- 
sured the inclination in all four bands (BVRl), which agreed 
with each other well within la errors (see Macri et al. 2000 for 
details), which wiU be used in the rest of the paper to cahbrate 
the TF relations and also to measure the value of Hq. 

In Appendix B, we summarize the inclination information 
compiled from previously published data for each galaxy in 
detail. In some cases, such as NGC 1365 or NGC 4725, the 
kinematical inchnation angle is quite different from the photo- 
metric inclination angles that we obtained. Although this could 
potentially lead to uncertainties in the corrected magnitudes and 
hnewidths, this should not introduce a bias in our analysis: for 
the distant cluster surveys, which will be used to measure the 
value of the Hq, the incUnations of galaxies were measured 
strictly from the imaging, either on photographic plates or on 
CCDs. Because our main objective is to create a calibration 
sample that is consistent with the cluster surveys, only the pho- 
tometric inclination angles (from Macri et al. 2000) will be 
considered for further analysis. 

2.3. Linewidths 
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Although linewidths have been measured and published for 
all the Key Project galaxies selected as calibrators for the Tully- 
Fisher relation, we have endeavored wherever possible to re- 
measure the various Unewidth definitions directly from the pro- 
files, in order to check for any anomalies. We are grateful to 
Martha Haynes for allowing us access to the profiles of the 
KP galaxies which were observed as part of the Giovanelli et 
al (1997a) database. In addition, several 21cm Une profiles 
were digitally scanned from the Uterature. We review here only 
briefly the linewidth database, as the details are explained in 
Macri et al. (2000). 

In order to avoid systematic effects when applying our 
Cepheid caUbration to various distant data sets, it is impor- 
tant that we use the applicable linewidth definition. The 20% 
linewidths, W20, are often used in the TP applications, which 
are the width in km s"^ of the HI profile measured at 20% of the 
peak of the HI flux, where for a 2-homed profile this is taken 
as the mean of each horn's peak flux. Another commonly used 
linewidth measurement is the 50% width, which is, for exam- 
ple, used in the 1-band data set of G97. The W50 is the width at 
50% of each of the horn maxima (or if the profile is Gaussian- 
shaped, then the width at 50% of the single peak). The potential 
advantage of measuring the linewidth at the 50% flux level is 
that it is less likely to be affected by noise, especially for low 
S/N profiles. However, the G97 Wso's are measured from an in- 
terpolation between the levels at 15% and 85% of the peak flux. 
For unusual profiles, this can result in underestimated profile 
widths. 

As all the KP galaxies are relatively nearby, most have very 
high S/N profiles. Many of them are also very large and overfill 
the telescope's beamsize. In these cases, to avoid systematically 
underestimating the profile widths, the observers have claimed 
to have been careful to ensure that a sufficient number of point- 
ings were obtained. 

Most galaxies have the classic steep-sided twin-homed pro- 
files of normal edge-on spirals. Unfortunately, there are ex- 
ceptions. Galaxies with different profiles are NGC 3031 
(which is distorted by a tidal interaction with M82=NGC 3034), 
NGC 3319 (which has a low-level pedestal), and NGC 3368 
(which has a non- symmetric two-homed profile). The 
widths are particularly sensitive to both low S/N profiles and 
pedestal features (which affect the positions of the 15% level) 
and unusual profiles (as the 85% level can be far from the end 
of the steep edge), the latter leading to systematically underes- 
timated widths. These cases are flagged by larger uncertainties. 

2.3.1. Corrections Applied to the Linewidths 

We apply the following corrections to the raw linewidth data. 
First, Unewidths were corrected for the inclination angle (i) and 
also for the redshift (z) effect following: 

W = Wismiy\l+zy^. (4) 

The inclination angle is derived from the eccentricity measured 
on the CCD frames via: 

_i j ib/af-ql 
V 1-^0 

where is the intrinsic minor-to-major axis ratio of edge- 
on spiral galaxies for which we adopt ^0 = 0.13 for T > 3 and 
qo = 0.20 for other types. No turbulence corrections are ap- 
pUed in this paper as the physical basis for this effect still re- 
mains ambiguous, except in estimating the intemal extinction 



corrections (in Section 3.2). Originally, one of the reasons for 
applying the turbulence correction was to straighten the TP re- 
lation, in addition to some theoretical preconceptions. Further- 
more, as we will discuss in later sections, the calibration galax- 
ies for our multi-wavelength TP relations mostly populate the 
region where the linewidths exceed ~300 kms"'. Since the 
turbulence correction is very small and constant for luminous 
galaxies (logl^ > 2.5), applying a turbulence correction would 
have no significant effect on our determination of Hq. 

2.4. Data Presentation 

In Table 2, we list all of the photometric, inclination and 
kinematical data used in calibrating the TP relation. The 
columns for Table 2 are: (1) name of the galaxy; (2)-(6) BVRI 
total and //-0.5 aperture magnitudes corrected for Galactic ex- 
tinction and internal extinction, as described in Section 2.2. For 
raw, uncorrected magnitudes, see Macri et al. (2000). The 
photometric data for NGC 224, NGC 598, NGC 2403 and 
NGC 3031 were taken directly from Pierce and Tully (1992); 
(7) logarithmic 20% linewidth corrected for inclination (us- 
ing the mean inclination) and redshift; (8) logarithmic 50% 
linewidth corrected for inclination and redshift; (9) photomet- 
ric I-band inclination; and (10) the mean photometric incU- 
nation which is used to calculate the internal extinction, and 
also to correct for inclination for the linewidths. We list here 
the I-band photometric inclination angles, since the strategy 
is to bring our photometry/Unewidth data to the same system 
as the distant cluster surveys. As will be described in detail 
later on, the I-band survey by G97 estimated the galaxy incli- 
nations from their 1-band images. We note, however, that in 
most cases, the BVRI photometric inclinations agree with each 
other extremely well, within la errors. 

3. CALIBRATION OF MULTI-WAVELENGTH TULLY-FISHER 

RELATIONS 

In this section, we derive TP cahbrations using two inde- 
pendent methods. The first is to obtain a consistent set of 
BVRIH^o s TP relations such that we can compare and examine 
the dispersions and the role of second parameters. This is done 
by using only the calibrator galaxies with directly-measured 
Cepheid distances. The 20% linewidths and the mean photo- 
metric inclination angles estimated from BVRI measurements 
will be used in this case. The second method is to derive TP cal- 
ibrations that can actually be applied to the distant cluster sur- 
veys to estimate Hq. This is done by using appropriate linewidth 
and incUnation angles, and also incorporating the cluster galaxy 
data to calculate the TF slopes. 

3.1. TF Calibration using Nearby Calibrator Galaxies Only 

First, we present BVRIH-o,i TP relations for 20% linewidth, 
using the mean of the BVRI photometric inclination angles (Col 
10 of Table 2), derived from 21 nearby galaxies with Cepheid 
distances only. While the optical photometric data, BVRI, are 
all total magnitudes, the H-band magnitude, H_o 5 is an aper- 
ture magnitude, extrapolated to the radius equivalent to half the 
total blue light of the galaxies. Slopes and zero points are de- 
termined using bivariate linear fits, minimizing errors in both 
logW andM'^: 



B^j. = -7.85(±0.7 l)(logW2o - 2.5) - 19.70(±0. 11) (6) 
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Table 2 

Photometric and Kinematical Data for lUlly-Fisher Calibrators 



Galaxy 










iic a 

-0.5 


logW^o' 






? 




mag 


mag 


mag 


mag 


mag 


kms"' 


kms"' 


deg 


deg 


N224'' 


-21.37 (19) 




-22.47 (18) 


-23.04(18) 


-23.83 (11) 


2.744(28) 


2.703(30) 


78(2) 


78(2) 


N598'' 


-18.57 (18) 




-19.39 (18) 


-19.80 (18) 


-20.31 (11) 


2.397(74) 


2.357(82) 


54(2) 


54(2) 


N925 


-19.59 (29) 


-19.95 (16) 


-20.24 (16) 


-20.59 (28) 


-21.16 (10) 


2.420(49) 


2.384(52) 


56(1) 


56(2) 


N1365 


-21.86 (37) 


-22.30 (12) 


-22.71 (17) 


-23.38 (12) 


-24.30(11) 


2.682(35) 


2.626(32) 


62(3) 


55(2) 


N1425 


-20.96 (17) 


-21.42(08) 


-21.82 (07) 


-22.38 (08) 


-23.07 (08) 


2.621(41) 


2.599(43) 


64(1) 


63(1) 


N2090 


-19.83(11) 


-20.20(10) 


-20.67 (10) 


-21.26(11) 


-21.95 (10) 


2.501(35) 


2.493(38) 


64(1) 


67(2) 


N24()3'' 


-19.11 (29) 




-19.93 (28) 


-20.32 (28) 


-21.14(25) 


2.480(59) 


2.457(62) 


58(2) 


58(2) 


N2541 


-18.72(18) 


-19.09 (13) 


-19.37 (10) 


-19.82 (13) 


-20.24 (10) 


2.370(49) 


2.361(52) 


58(2) 


62(1) 


N3031" 


-20.64(17) 




-21.87 (17) 


-22.43 (17) 


-23.54 (10) 


2.719(34) 


2.654(39) 


58(2) 


58(2) 


N3198 


-20.28 (08) 


-20.75 (08) 


-21.07 (07) 


-21.55 (08) 


-22.22 (08) 


2.531(32) 


2.506(34) 


71(2) 


68(2) 


N3319 


-19.34(14) 


-19.69 (14) 


-19.94 (13) 


-20.32 (14) 


-20.44 (13) 


2.4Uj(4iS; 


2.j64pl) 


60(1) 


58(1) 


N3351 


-19.78 (11) 


-20.56 (09) 


-21.00(09) 


-21.61 (09) 


-22.56 (10) 


2.586(47) 


2.577(42) 


45(2) 


45(2) 


N3368 


-20.47 (14) 


-21.26(12) 


-21.73 (11) 


-22.29(11) 


-23.41 (11) 


2.674(36) 


2.631(34) 


51(1) 


49(1) 


N3621 


-19.48 (12) 


-20.00(12) 


-20.50 (12) 


-20.99 (12) 


-21.84 (12) 


2.499(35) 


2.471(38) 


64(1) 


66(1) 


N3627 


-21.10(18) 


-21.73(18) 


-22.12 (18) 


-22.63 (18) 


-23.50 (18) 


2.626(26) 


2.694(35) 


65(1) 


65(1) 


N4414 


-20.88 (13) 


-21.57(11) 


-21.98(11) 


-22.59(11) 


-23.65 (11) 


2.743(39) 


2.562(39) 


46(2) 


46(2) 


N4535 


-20.80 (10) 


-21.38 (09) 


-21.61 (08) 


-22.27 (08) 


-22.72 (09) 


2.586(38) 


2.538(32) 


49(1) 


51(1) 


N4536 


-20.49(12) 


-21.05 (10) 


-21.39 (09) 


-21.95(13) 


-22.79 (09) 


2.562(30) 


2.695(31) 


68(1) 


69(1) 


N4548 


-20.36 (24) 


-21.04(23) 


-21.47 (23) 


-22.13(23) 




2.617(46) 


2.729(25) 


38(1) 


39(1) 


N4725 


-21.27 (10) 


-21.85 (09) 


-22.20 (09) 


-22.74 (09) 


-23.65 (10) 


2.671(26) 


2.564(30) 


51(1) 


62(1) 


N7331 


-21.58 (12) 


-22.30 (11) 


-22.73 (11) 


-23.29(11) 


-24.66 (11) 


2.746(21) 


2.576(45) 


68(2) 


69(2) 


"The fomat of magnitudes is sucli tliat -21.65(13) is equivalent to -21.65 ± 0.13. 










""Tlie format of log W is such that 2.682(35) is equivalent to 2.682 ± 0.035. 













'^The inclination errors are shown in brackets. 



•^Magnitudes from Pierce and TXilly 1992 



= -8.81(±0.82)(logW2o-2.5)-20.27(±0.12) (7) 

=-8.68(±0.71)(logW2o-2.5)-20.60(±0.11) (8) 

/^ = -9.21(±0.75)(logW2'i,-2.5)-21.09(±0.12) (9) 

//f()5 =-11.03(±0.86)(logW2o-2.5)-21.74(±0.14) (10) 

In Figure 1, the above TF relations are plotted. The observed 
dispersions of above relations are a = 0.45,0.37,0.35.0.37 and 
0.36 for BVRI and H_o.5 respectively, which are combinations 
of observational errors and intrinsic scatter of the TuUy-Fisher 
relations themselves. The latter will be discussed in Section 
3.4. Errors in logW and absolute magnitudes were determined 
as follows. The linewidth error is a combination of two sources: 
the observational error in the Unewidth measurement and the 
uncertainty in the inclination angle. For absolute magnitudes, 
their uncertainties come from a few sources including the un- 
certainty in the Cepheid distance modulus (random only, typ- 
ically ~0.10 mag), observational photometric error (typically 
~0.10 mag), the error in internal extinction law which is asso- 
ciated with the error in the inclination angle, and the uncertainty 
in the Galactic extinction (^0.02 mag in 1-band). 

We have also calculated the slope and zero point of each TF 
relation using the direct and inverse fits for each wavelength 
as a consistency check, in which the fits were determined by 
minimizing the errors in magnitudes and linewidth only respec- 
tively. Schechter (1980) suggested that the average value of 
\ogW at a constant magnitude does not depend on whether the 
sample is distance- or magnitude-Umited. That is, the inverse 
TF relation is insensitive to the incompleteness bias. Based on 
this, some authors have preferred to use the inverse TF rela- 
tion to derive cluster distances (e.g. Tully 1999). Although the 
slopes do change sUghtly as expected depending on what fit is 
used, the zero point remains well within \a. As wiU be shown 



in Section 5.2.1, the TF slope uncertainty has a very little effect 
on the Ho value. In this paper, we will use the bivariate fits only. 
Several articles have investigated the possibility of a TF relation 
being a quadratic relation (cf: Mould, Han and Bothun 1989). 
We do not use quadratic TF relations in the main analyses of 
this paper, however, interested readers are referred to Appendix 
C for further discussion. 

3.2. TF Calibration using Both the Calibrators and Cluster 

Galaxies 

We can also calibrate TF relations by combining the local 
calibrator galaxies and cluster survey galaxies. These indepen- 
dent samples should supplement each others' advantages and 
disadvantages. The set of equations (6-10) will provide a con- 
sistent set of TF calibrations which can be used to investigate 
several fundamental questions pertaining to the TF correlation 
itself, which will be discussed in the following sections. How- 
ever, in order to fulfill the original goal of this project, to mea- 
sure the value of Hq, we must in addition derive TF equations 
that are exactly on the same photometric-linewidth system as 
the distant cluster surveys. Thus, a second set of the / TF rela- 
tions was derived using slightly different methods. 

The G97 I-band survey used I-band photometric inclina- 
tions. In determining the zero points and slopes, we take ad- 
vantage of the all-sky survey data that are available. Calibrat- 
ing the TF relations using only local galaxies with Cepheid dis- 
tances poses two problems. This sample, which consists of only 
21 galaxies, is not statistically large enough to constrain the 
slope very well, and errors in the individual galaxy distances 
contribute significantly to the final dispersion. To minimize the 
first problem, we iteratively determine the slope and zero point 
using both the distant cluster sample and the local calibrators 
as follows: (1) the slope and zero point are calculated using 
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the calibrator sample only. (2) Then distances to distant cluster 
galaxies (e.g. the G97 1-band sample) are determined using 
the caUbration from (1), such that all the clusters are put on the 
same (log W- 2.5) -M/ plane. Here, the cluster sample galaxies 
are corrected for the cluster incompleteness bias (which will be 
discussed in detail in Section 5.4). (3) The slope is determined 
using this combined cluster sample. (4) The zero point is then 
estimated by minimizing the dispersion in the local calibrating 
sample using the slope from (3). Steps (2) through (4) are re- 
peated until the slope and zero point converge. Convergence is 
usually achieved within three or four iterations. We obtain the 
following I-band TF relation: 

7J =-10.00(±0.08)(logW5'o- 2.5) -21.32 (11) 

For reference, the I-band-50% TF relation derived us- 
ing only the calibrators (equivalent to Equation 9) is: Ij = 
-9.87(logW5Q- 2.5) -21.30, which is virtually the same as 
Equation 1 1 . Finally, note that the internal extinction correc- 
tion of T98 is expressed in terms of 20% Unewidths. The 
G97 database, on the other hand, lists 50% Unewidths. We 
have made a compilation of galaxies with both published 20% 
and 50% Unewidths by examining the G97 survey and an H- 
band survey by Aaronson, Huchra, Mould and collaborators (cf: 
Aaronson et al. 1982). The mean ratio for 73 such galaxies is 
W20 = 1 . l(±0.03)W5o. When estimating the internal extinction, 
we calculate the corresponding 20% linewidth using this ratio 
for all the cluster galaxy data. In addition, the T98 internal 
extinction expressions are given as a function of turbulence- 
corrected linewidth. Thus, only in estimating the internal ex- 
tinction, we apply this turbulence correction to the linewidths 
following Tully & Fouque (1985). 

3.3. Uncertainties in the Tully-Fisher Parameters 

To better understand the errors in the derived TF relations, 
we generated a sample of 5000 random points which followed 
a Tully-Fisher relation of dispersion 0.43 mag. From this sam- 
ple, 500 random subsamples of galaxies were drawn, and the 
zero point and slope were calculated for each trial run. We ex- 
amined three cases: N=25, 500 and 2000. For A' = 25 trial runs, 
the range of slope distribution is -^5 times larger than other two 
larger samples, suggesting that the sample solely consisting of 
the ^20 Local calibrators could give a slope estimate that is 
significantly different from the true value. The cluster database 
used to estimate the TF slope was comprised of more than 500 
galaxies (for the I-band survey). Even for this larger sample, 
the slope distribution has a \a RMS scatter of 0.09 mag. 

The distributions of zero points obtained for the same set of 
simulations were also examined. Again as expected, the small- 
est sample, in which 25 galaxies were picked during each trial 
run, has the largest standard deviation in the zero point distri- 
bution of 0.13 mag, which in fact agrees with the actual zero 
point error calculated in Equations 6-10. 

We have detemuned the slope using a sample of ^^500 galax- 
ies and then the zero point using a sample of ~20 galaxies. 
Thus, we formally adopt the uncertainties for the slope and 
zero point of 0.09 and 0.13 mag respectively. In Section 5.5, 
the sensitivity of the value of Hq to the slope uncertainty of the 
TF calibration will be investigated further. 

3.4. Dispersion 

The dispersion of the TF relation has three primary origins: 
(1) intrinsic to the relation itself, due to departures from the 



"perfect" physical correlations, such as Freeman's exponential 
disk law, or variations in mass-to-light ratio, (2) observational 
errors in magnitude and Unewidths, and (3) for the case of a 
cluster TF relation, the physical depth of the cluster itself or, 
for the case of the calibrators, the uncertainties in the distance 
estimates of the individual galaxies. 

The subject of the TF dispersion has long been of interest, as 
its value can pose strong constraints on scenarios for galaxy for- 
mation and evolution. For example, Eisenstein and Loeb (1996) 
suggested that a dispersion exceeding 0.3 mag is expected for 
most galaxy formation histories. Silk (1996) speculated that a 
tight TF relation could suggest an unexpectedly late epoch of 
galaxy formation. It is also important to understand the intrin- 
sic scatter of the TF relation, because it is directly related to the 
bias in the distance ladder itself. 

Surprisingly, estimates of dispersion have been very ambigu- 
ous, ranging from ^ 0.1 mag (e.g. Bernstein et al. 1994, Freed- 
man 1990), up to ~ 0.7 mag (e.g. Kraan-Korteweg et al. 1988, 
Sandage 1994, Tammann 1999), with many estimates around 
0.2-0.3 mag (e.g. Pierce & Tully 1988, WiUick 1991, Courteau 
1992, Bothun & Mould 1987, Schommer et al. 1993). More 
recent studies have shown furthermore that the TF dispersion 
is a more complex function of galaxy types and sizes. Using 
his LCO/Palomar Survey, WiUick (1999) found that the disper- 
sion is as much as ~0.75 mag for low-luminosity, low-surface 
brightness galaxies, while for the higher luminosity, higher sur- 
face brightness galaxies, it is < 0.35 mag. GiovanelU (1997) 
reported that the dispersion of the TF relation is in fact a func- 
tion of linewidth; it decreases from fast-rotators (cr ^ 0.2 mag 
for logW ~ 2.6) to slow-rotators (> 0.3 mag for logW < 2.2). 
Readers are referred to GiovaneUi (1997) and WiUick (1999) for 
an extensive discussion on this subject. Here, we have a unique 
set of data: the local calibrator data provide a sample in which 
the distance to each galaxy is known with 10-15% accuracy un- 
like the case of the cluster TF relation, where the cluster depth 
can be an additional unknown parameter. Freedman (1990) also 
used galaxies with accurately measured Cepheid distances, but 
this sample now has 21 galaxies, instead of five used by Freed- 
man (1990). Dispersions forBVRIH Tully-Fisher relations ob- 
tained from our Cepheid calibrator galaxies were quoted in Sec- 
tion 3.1. These are, however, observed dispersions and are a 
combination of observational errors and the intrinsic scatter of 
the TF relation itself, added in quadrature. We determined the 
intrinsic scatter of the TF relation for each wavelength, by esti- 
mating the observational error for each calibrator galaxy fol- 
lowing Equation 10 of GiovanelU et al. (1997b). We have 
also assumed that the dispersion is a constant over the linewidth 
range 2.4 < logl^ < 2.8 spanned by the calibrators. A typical 
observational error for galaxies is ^ 0.30 ± 0.05 mag, (whose 
main sources are photometric, distance modulus uncertainty 
and the linewidth measurement uncertainty multiplied by the 
TF relation slope,) indicating that the intrinsic dispersions are 
0.25,0.22,0.25,0.20 and 0.19 mag for B, V, R, I, and H_o.5 re- 
spectively. We have tested to see if this dispersion is a function 
of some observational parameter, such as the color. In Figure 
2, we plot TF residuals as a function of color, which show no 
dependence, suggesting that color is certainly not a candidate 
for the second parameter in TF relations. 

Another question to be addressed is whether a dispersion of 
0.20-0.25 mag is a value that could have been measured acci- 
dentally. In the previous section, the slope and zero point of 
the TF relation were examined using a randomly created sam- 
ple of 5000 galaxies. A sinular exercise is carried out here. 
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except a sample of 5000 points were generated assuming that 
the dispersion of the TF relation is 0.7 mag. A random subsam- 
ple of galaxies was again withdrawn 500 times. The dispersion 
of 0.2 mag is highly unlikely, with a probability of practically 
zero. The B- and V-band TF relation dispersions were slightly 
higher than that of other wavelengths, measured at 0.25 mag. 
The chance of 'accidentally' observing even this dispersion is 
remote, less than 1%, if the true dispersion is 0.7 mag. 

Furthermore, when the cluster data, consisting of few hun- 
dred galaxies, are put on the same TF plane, one finds that its 
dispersion is less than 0.3 mag. This confirms that the intrinsic 
dispersion of TF relation is relatively small. 

3 .5 . Color Dependence of Tully-Fisher Relations 

If the internal extinction correction were perfect, there should 
be no dependence of the corrected color of the TF galaxies on 
inclination. In Figure 3 (top), we examine the (B-I) color of 
the local calibrators as a function of inclination. The solid cir- 
cles represent the data corrected using Equation I, while the 
open circles indicate the uncorrected data. There is a slight dip 
in this color-inclination relation around the inclination angle 
of 50°, however, there is no overall trend that is observable. 
Also shown in the middle and bottom plots of Figure 3 are the 
{B - 1) vs inclination relations for the same set of Local cali- 
brator galaxies, but here, Han's (Equation Al in Appendix A) 
and Giovanelli et al.'s (Equations A2 and A3) extinction cor- 
rections, respectively, were used. Although for both Han's and 
Giovanelli et al.'s corrections there may be a slight decline in 
the relation for the large-inclination region, this apparent trend 
could be partly due to small-number statistics. One significant 
question is whether the value of Ho is affected by the uncer- 
tainty arising from the internal extinction corrections. This wiU 
be discussed in detail in Section 5.7. 

3.6. Barred Galaxies in the Tully-Fisher Relations 

Among the 21 calibrator galaxies in our Tully-Fisher sample, 
seven are barred galaxies. In this section, we explore whether 
these barred galaxies are in fact suitable TF caUbrators or not. 
On the left side of Figure 4, we have re-plotted B, I and H 
band TF relations using the 20% linewidth. Seven barred galax- 
ies, NGC 925, NGC 1365, NGC 3319, NGC 3351, NGC 3627, 
NGC 4535, and NGC 4725 are shown with open circles. The 
solid line in each plot represents the TF fit to all the galaxies. 
A striking feature is that the majority of them lie above the TF 
relations, suggesting that they are systematically brighter for a 
given rotational velocity. This is true in all three wavelengths. 
(Although not shown here, the same applies to the V- and R- 
band TF relations). In Table BI, a summary of published pho- 
tometric and kinematical inclination angles are listed. We use 
the kinematical inclination angles for the seven barred galaxies 
and plot the TF relation on the right side. The soUd line in each 
plot shows the fit through all the galaxies, using the kinematical 
inclination angles for the barred galaxies. The dotted line is the 
fit through the TF relation using the photometric inclination an- 
gles for all galaxies (same as the solid line in the left plot). The 
change in zero point in the I-band, for example, is 0.04 mag. 
Also, for the I-band, the observed dispersion decreased con- 
siderably, from 0.37 mag down to 0.32 mag. For the B-band 
relation, it changed significantly as well from 0.45 to 0.41 mag. 
The dispersion remains unchanged when kinematical inclina- 
tion angles are substituted for all other galaxies, indicating that 
the photometric determination of inclination angles for barred 



galaxies are systematically different from the kinematical esti- 
mates. 

Two possible explanations for the significant difference in 
the measured photometric and kinematical inclination angles of 
barred galaxies. The first is that the star formation rate is higher 
in the barred galaxies, although observations suggest that the 
effect of a bar on the disk star formation rate is unimportant 
(Kennicutt 1999). In Figure 5, colors of the calibrator galaxies 
are shown as a function of linewidth. The barred galaxies are 
indicated by open circles. In the plots on the left side, the pho- 
tometric inclination angles were used to derive the corrected 
magnitudes, while on the right side, the magnitudes were cor- 
rected using the kinematical inclination angles for seven barred 
galaxies. Although most barred galaxies lie on the blue edge 
of this correlation if the photometric inclination angles are used 
(with one exception, NGC 3351, which is located at the red 
edge), when the kinematical inclination angles are substituted 
for barred galaxies, there are no longer any obvious differences 
between the distributions of the barred and non-barred galax- 
ies. If an additional star formation is actually triggered by the 
presence of a bar, then one would expect the deviation of barred 
galaxies from the TF relation would be higher for shorter wave- 
length. We do not see such evidence, although our sample is 
rather small. We find the deviations from the mean TF fit for 
barred galaxies are the same for all five wavelengths. The sec- 
ond explanation, which is more plausible, is that the incUnation 
angles for barred galaxies are over-estimated if the photometric 
inclination angles are used. This would make the barred galax- 
ies appear brighter for a given linewidth, and would explain that 
adopting the kinematical incUnation angles would shift these 
galaxies closer to the TF relation. 

As stated in Section 2.2.2, in the following analysis in this 
paper, we will not use the kinematical inclination angles, be- 
cause those for the I-band cluster survey galaxies were de- 
rived strictly from the photometric observations. In Figure 6, 
we show a TF relation for cluster galaxies. Those classified 
as barred are shown by solid circles. There is no obvious sys- 
tematic offset between the barred and non-barred samples here, 
although one would expect that both the morphology identifi- 
cation and inclination estimates might be more difficult for dis- 
tant galaxies. We have also examined the distribution of the 
ratio of number of barred to non-barred galaxies as a function 
of linewidth or inclination, but observed no significant differ- 
ence between the cluster and caUbrator samples. 

4. APPLICATION OF THE TULLY-FISHER RELATION TO DISTANT 
CLUSTER SAMPLES AND THE VALUE OF THE HUBBLE 
CONSTANT 

There have been numerous Tully-Fisher surveys of clusters 
of galaxies in various wavelengths. Most recently, Giovanelli 
and collaborators presented a survey of I-band photometry and 
radio hnewdiths of ~2000 spiral galaxies in clusters out to 
~ 10,000 km s~^ Aaronson, Huchra, Mould and collabora- 
tors presented H-band aperture magnitudes and 20% linewidths 
of galaxies in clusters out to ^15,000 km s"' (Aaronson et al. 
1980, 1982, 1986 : hereafter AHM). Bothun et al. (1985) pre- 
sented a catalog of UBVR multi- aperture photometry, H-band 
photometry and 21cm HI observations for several hundred spi- 
ral galaxies in 10 clusters. Han and Mould (1990) published 
I-band magnitudes for galaxies in the AHM survey. In this pa- 
per, we focus on the I-band survey by Giovanelli et al. (1997a), 
as it is the most complete one to date, and presents total I-band 
magnitudes measured from CCD observations, unlike in other 
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surveys where aperture magnitudes were used. 

4. 1 . I-band Cluster Distances and the Hubble Constant 

The absolute calibration of the BVRIH-0.5 TuUy-Fisher rela- 
tions was presented in the previous section. The I-band rela- 
tion was also re-calibrated for specific application to the G97 
database. The second I-band calibration can now be applied 
to the distant cluster surveys, which extend out to --^ 10,000 km 
s~\ to determine Hq. 

For a detailed description of the I-band database, readers are 
referred to GiovanelU et al. (1997a) who explain in detail the se- 
lection criteria and the cluster galaxy membership. Briefly, the 
G97 survey is comprised of clusters and groups of galaxies out 
to 9000 km s~^ . The database consists of 50% linewidths, to- 
tal /-band magnitudes which were extrapolated to infinity (the 
largest isophotal elUpses measured were usually ~24.5 mag 
arcsec"^), and photometric inclination angles measured from 
I-band images. 

In Table 3, we summarize the clusters of galaxies that are 
used in the following analyses to measure the Hubble constant. 
Also Usted in this Table are three velocities for each cluster: ve- 
locities corrected for the Local Group frame, cosmic microwave 
background frame (CMB), and corrected for the flow model 
which will be described in Section 5.8. 

To be included in the TF analyses, the galaxies had to meet 
all of the following criteria: 

(1) They do not deviate by more than twice the dispersion 
from the TF relation. Three galaxies were excluded due to this: 
11830 in the Fornax cluster, N2535 in Cancer, and N3832 in 
A1367. 

(2) Galaxies must not be nearly face-on (; < 40°) because 
for these objects, the inclination uncertainty and subsequently 
the linewidth uncertainty is too large. Also, galaxies should not 
be close to edge-on (i > 80°) since the extinction correction for 
these become unreasonably large. 

(3) The linewidth distribution of galaxies in the clusters 
should be similar to that of the calibrators. In Figure 7, some 
examples of clusters are shown. In the same figure, the cal- 
ibrator data are over-plotted on each cluster TF relation with 
solid circles. The local calibrators and the cluster samples do 
not all cover the same rotational velocity range. For some clus- 
ters, such as Eridanus, Cancer, or Pegasus groups, the cluster 
sample extends to much smaller linewidths. In order to make a 
consistent sampling between the two sets of datasets, we apply 
a lower logW cutoff of 2.35. 

(4) The internal extinction correction should not exceed 0.6 
mag. In Figure 8, the distributions of internal extinction values 
for the cluster galaxies and calibrators are obtained using three 
separate extinction corrections, that of T98, Han and Giovanelli 
et al. (see Appendix A for the latter two corrections). For the 
T98 and Giovanelli et al. corrections, there is a significant dif- 
ference between the two samples: the cluster galaxies seem to 
have a long upper-end tail in their distribution of extinction val- 
ues. To make a consistent sampling between the calibrators and 
clusters, we apply an upper cutoff of A,„, / = 0.6 mag. It is also 
a matter of common sense to exclude galaxies more than half 
of whose Ught is obscured (Ai„, > 0.75 mag). 

The final /-band sample, referred to as the G97 sample, thus 
consists of 276 galaxies. 

We measured the TF distances to all the clusters with 5 mem- 
ber galaxies or more by applying Equation 1 1 to the I-band sur- 
vey. The results are shown in Tables 4 in which the cluster or 
group distances, D, and V/D values (corresponding to the ve- 



locities with respect to the Local Group, CMB and model) are 
listed for 1-band clusters. The uncertainty in the Hq values is 
roughly a combination of observational errors, errors in the TF 
zero points and slopes, and the dispersion of the TF relation. 
This will be discussed in more detail in Section 6. 

The top plot in Figure 9 shows the Hubble diagram for the 
I-band clusters. Also in the bottom of the same figure, we plot 
the Hubble constant for each cluster as a function of velocity. 
In order to measure the values of Hq, we use only those clus- 
ters satisfying (1) Vcmb > 2000 km s~', and (2) velocities with 
respect to the CMB and the flow field model are within 10% of 
each other. This leaves us with a set of 15 clusters and groups 
of galaxies. Also, it will be shown in Section 5.1.2 that there 
are no selection biases for galaxies with Vcmb > 2000 km s~', 
thus using those clusters with Vcmb > 2000 km s"' should be a 
reasonable choice. Taking an average of 15 clusters, we obtain 
Ho = 73 ± 2 (random) ± 9 (systematic) using the I-band sur- 
vey. The random and systematic errors are discussed in detail 
in Section 6. 

4.2. Comparison with Other Surveys 

As mentioned above, there are other surveys of clusters of 
galaxies that reach velocities ^10,000 km s"', which are suit- 
able for the determination of Hq. The G97 survey was used as 
the primary survey in this paper because of its completeness and 
quality. In this section, we compare our I-band Hq results with 
the distances and values of Hq obtained using other surveys. 

4.2.1. H-band 

The H-band cluster surveys come from two sources: Aaron- 
son et al. (1982) and Aaronson et al. (1986). We refer to 
the combined sample as AHM. In these surveys, rotational ve- 
locities were measured at the 20% level, and //-band aper- 
ture magnitudes were corrected to an isophotal aperture of 
logA/Di = -0.5 by adopting the system of the RC2 catalog (de 
Vaucouleurs, de Vaucouleurs, and Corwin 1976). These magni- 
tudes are referred to as //-0.5 (see Aaronson, Mould and Huchra 
1980 for details of how the interpolation was made). Also, in 
the original database of Aaronson et al. (1980), they included 
a 3° additive term for inclination angles. We are not adopting 
this convention, because we use the axis ratios given by RC3 
catalog. 

Cluster membership in the H-band survey was not as well 
defined as for the I-band survey, where G97 thoroughly in- 
spected each cluster to determine which galaxies are likely 
members of the cluster. We examined the distribution of galax- 
ies in each cluster in both spatial and velocity space to identify 
any outiiers. Several clusters, including A1367, Z74-23, Her- 
cules and Pegasus, were split into two or three subgroups. 

The same four criteria listed above for the selection of the 
I-band survey galaxies are applied to the H-band AHM sur- 
veys. Again, to make the cluster database consistent with the 
calibrators, a lower cutoff of log W > 2.35 is applied. The final 
H-band sample thus consists of 163 galaxies in 26 clusters. 

Following the method described in Section 3.2, the H-band 
TF relation was derived using both the calibrators and the clus- 
ter data. The resulting value of Hq is 67 ± 3 (random) ±10 
(systematic). The 10% discrepancy between the H-band and I- 
band values of Ho is statistically significant, greater than the 2(j 
level. Tormen & Burstein (1995) published a list of isophotal H 
magnitudes for galaxies from the AHM survey, corrected using 
RC3 parameters. Substituting this set of data for both the cal- 
ibrators and cluster galaxies yields the same value of Hq; use 
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Table 3 

Clusters of Galaxies in the I-band Survey 



Cluster/ 


R.A.' 


DEC" 




y 


^ Model 


N 


Group 


hh mm ss 


dd mm ss 


km s 


km s 


km s 




N383 


01 04 30 


+32 12 00 


5425 


4924 


5086 


23 


N507 


012000 


+33 04 00 


5350 


5869 


5016 


14 


A262 


01 49 50 


+35 54 40 


5170 


4730 


4852 


32 


A400 


02 55 00 


+05 50 00 


7258 


7016 


6983 


26 


Eridanus 


03 3000 


-21 30 00 


1641 


1607 


1627 


35 


Fornax 


03 36 34 


-35 36 42 


1328 


1380 


1372 


44 


Cancer 


08 17 30 


+21 14 00 


4652 


4982 


4942 


29 


Andia 


10 27 45 


-35 04 1 1 


2505 


3106 


2821 


26 


Hydra 


10 34 28 


-27 16 26 


3454 


4061 


3881 


27 


N3557 


1 1 07 35 


-37 16 00 


2702 


3294 


2957 


12 


A 1367 


11 41 54 


+20 07 00 


6316 


6709 


6845 


38 


Ursa Major 


1 1 54 00 


+48 53 00 


957 


1088 


1088 


30 


Cen30 


12 46 06 


-41 02 00 


2756 


3272 


4445 


39 


Cen45 


12 47 57 


-40 22 23 


4304 


4820 


4408 


9 


Coma 


12 57 24 


+28 15 00 


6889 


7143 


7392 


41 


ESO508 


13 09 54 


-23 08 54 


2662 


3149 


2896 


17 


A3574 


13 46 06 


-3009 00 


4308 


4749 


4617 


20 


Pavo 


20 13 00 


-71 00 00 


3919 


4027 


4220 


10 


MDL59 


22 00 18 


-32 14 00 


2621 


2304 


2664 


26 


Pegasus 


23 17 43 


+07 55 57 


4109 


3545 


3874 


19 


A2634 


23 35 55 


+26 44 19 


9516 


8930 


9142 


27 



"Positions given in Epoch 1950.0. 



Table 4 

I-band TUlly-Fisher Distances to Clusters 



Cluster 


N 


m—M 


D 




^cmb/o 


Vodel/O 






(mag) 


(Mpc) 


(km/s/Mpc) 


(km/s/Mpc) 


(km/s/Mpc) 


A1367 


28 


34.84 


93.0 


67.9 (11.0) 


72.2 (11.7) 


73.6 (11.9) 


A2197 


3 


35.56 


129.3 


71.8(11.4) 


70.5(11.1) 


73.9 (11.7) 


A262 


22 


34.18 


68.7 


75.3 (12.2) 


68.9 (11.1) 


74.2(12.0) 


A2634 


20 


35.36 


118.0 


80.6 (12.6) 


75.7 (11.8) 


77.5 (12.1) 


A3574 


14 


34.03 


64.1 


67.2 (10.5) 


74.1 (11.6) 


72.0(11.2) 


A400 


18 


34.81 


91.8 


79.1 (12.3) 


76.5(11.9) 


76.1 (11.8) 


Antlia 


16 


33.30 


45.6 


54.9 ( 8.8) 


68.1 (10.9) 


61.8(9.9) 


Cancer 


17 


34.40 


75.8 


61.4(9.8) 


65.8 (10.5) 


65.2(10.4) 


Cen30 


22 


33.25 


44.7 


61.7(10.2) 


73.2(12.1) 


99.5 (16.4) 


Cen45 


8 


34.24 


70.5 


61.1 (10.0) 


68.4(11.2) 


62.5 (10.3) 


Coma 


28 


34.74 


88.6 


77.7 (12.2) 


80.6(12.6) 


83.4(13.0) 


Eridanus 


14 


31.66 


21.5 


76.4(12.4) 


74.8 (12.1) 


75.8(12.3) 


ESO508 


9 


33.07 


41.1 


64.8(10.3) 


76.7 (12.1) 


70.5 (11.2) 


Fornax 


14 


30.93 


15.3 


86.6(14.0) 


90.0(14.5) 


89.5 (14.4) 


Hydra 


17 


33.90 


60.2 


57.4 ( 8.9) 


67.4(10.5) 


64.5 (10.0) 


MDL59 


10 


32.56 


32.5 


80.7 (12.6) 


71.0(11.1) 


82.1 (12.8) 


N3557 


9 


33.02 


40.2 


67.2(11.1) 


81.9(13.5) 


73.5 (12.1) 


N383 


15 


34.19 


68.9 


78.7 (12.3) 


71.5(11.2) 


77.3 (12.1) 


N507 


7 


33.86 


59.1 


90.6(14.0) 


82.4(12.8) 


89.0(13.8) 


Pavo 


5 


32.71 


34.8 


112.6(17.4) 


115.7(17.9) 


121.3 (18.7) 


Pavo2 


13 


33.62 


52.9 


81.3(13.0) 


83.2(13.4) 


87.9 (14.1) 


Pegasus 


14 


33.73 


55.8 


73.6(11.5) 


63.5 ( 9.9) 


69.4 (10.9) 


Ursa-Major 


16 


31.58 


20.7 


46.3(7.1) 


52.6(8.1) 


52.6(8.1) 
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of these updated magnitudes does not reconcile the disagree- 
ment between 1-band and H-band Ho values. We defer further 
discussion of this discrepancy to Section 4.2.3. 

4.2.2. B- and V-band Surveys 

Bothun et al. (1985: hereafter B85) presented a catalog of 
21cm HI observations, UBVR multi-aperture photometry, and 
H-band photometry of several hundred galaxies in 10 clusters. 
The total UBVR magnitudes were extracted from multiple- 
aperture measurements by applying the aperture corrections es- 
timated from the mean growth curve derived from the RC2. The 
average aperture correction amounted to 0.2-0.3 mag. Of these 
10 clusters, we selected five clusters, A1367, A400, A539, 
Coma and Z74-23, which consisted of at least five galaxies 
with good quality radio and optical data. Using 57 galaxies 
with logWio > 2.4 in five clusters, we obtain Ho = 80 ± 8 km 
s"' Mpc"^ (random) for B-band, and Ho = 73 ± 6 km s"^ Mpc"^ 
(random) for V-band. However, we note that the Ho value for 
A 539 is almost 3(t larger than those of other four clusters, in 
both B-band and V-band estimates. Excluding this cluster, we 
obtain Hq = 72 ± 2 km s"^ Mpc"' for B-band and 68 ± 2 km 
s~' Mpc~' for V-band. These values lie between the I-band 
and H-band results. We would put less weight on the B-band 
and V-band magnitudes as they were extrapolated from aper- 
ture magnitudes, and our B- and V-band calibrator magnitudes 
were not derived in the exactly the same manner. 

4.2.3. Discussion 

We have calculated Hq using clusters of galaxies in four 
wavelengths: B, V, I and H_o.5. The derived values of Hq vary 
from 67 up to 73 km s~' Mpc~'. In this section, we focus on 
the difference between the 1-band and H-band results which 
are based on more reliable databases, and are statistically better 
as the number of cluster galaxies used to derive the value of Ho 
is more than twice as large as those in the B- and V-band sur- 
veys. It appears that the H-band TF relation yields consistently 
larger distances. We examine the databases carefully with other 
available published data in order to understand the root of the 
Ho discrepancy. 

Linewidth: The I-band analysis is based on a different 
linewidth set (W50) than the H-band analysis (W2()), so we first 
tested whether the discrepancy in distances was caused by a 
subtle error in one of the linewidth sets. Because TF rela- 
tions are steep, it is important to measure linewidths accurately. 
Bothun et al. (1985) presented both 50% and 20% linewidths, 
and the 50% widths agree well with the published values of the 
G97 sample: < W^Bothun/^G97 >= 1-000 ±0.003. As noted 
in Section 2, we re-measured 20% and 50% linewidths for the 
local calibrators, and found that the W20/W50 ratios are consis- 
tent with those in the cluster samples, W2Q /W50 = 1 .07 ± 0.03 
and 1.10 ±0.05 respectively. Furthermore, a list of galax- 
ies that are in both G97 and AHM samples was compiled. 
The mean W20/W50 ratio of galaxies in the overlapping sam- 
ple (W20/W50 = 1.10±0.08) again agrees well with the mean 
ratio of the local calibrators. We have also compiled a set of 
-^70 galaxies which had both I-band and H-band magnitudes, 
and also 20% and 50% linewidths measured, and confirmed that 
the derived values of Ho are insensitive to whether 20% or 50% 
linewidths are used. Han & Mould (1992) presented an I-band 
TF survey that used 20% linewidths. As an additional check we 
repeated the entire analysis using the Han & Mould database. 
We obtained distances that are consistent with those derived 
above using the G97 sample. This leads us next to the various 



tests to confirm the photometric quality of the cluster data and 
calibrators. 

Photometry: We first tested the zero point of the G97 I-band 
magnitudes by comparing them to those published by Han & 
Mould (1992). There were 59 galaxies that overlapped in the 
two surveys: they agree within 0.01 mag in the mean. On the 
other hand, Peletier and Willner (1993) published //-0.5 mag- 
nitudes for Ursa Major cluster galaxies. We found that these 
also agreed within la error with the AHM magnitudes. The 
R-I colors of the TF calibrator galaxies (whose magnitudes 
are given by Macri et al. 2000) agree well with those of normal 
spiral galaxies presented in de Jong (1998). 

Maeder et al. (priv. comm.) show that the concentric aper- 
ture photometry of Aaronson et al is verified by recent 2MASS 
survey photometry. But the consistency between calibrators' 
and clusters' isophotal diameters, on which H_o.5 magnitudes 
are based, is more problematic. Diameters on the RC2 and RC3 
systems differ by 0. 1 dex (Tormen & Burstein 1995), and, while 
this apparently applies both to caUbrators and cluster galaxies, 
the absence of systematic differences of the same order (equiv- 
alent to 0.2 mag in H_o.5) cannot be guaranteed without further 
work. 

Color Distribution of Calibrators and Cluster Galaxies: 

The results above suggest that the difference in 1-band and H- 
band distances must reflect a systematic difference in measured 
I-H-0,5 colors between the calibrator and cluster samples. In 
the top panel of Figure 10, the I-band and H_o.5-band TF dis- 
tance moduli are compared. In the middle plot, the /-//_o.5 
color of the same galaxies is shown as a function of linewidth. 
The calibrators are also shown by solid circles. The most strik- 
ing feature is that calibrators are redder compared to the cluster 
sample galaxies (see below), which is the most likely reason 
the H- and I-band TF distances disagree at a 10% level (top 
figure). Environmental effects could possibly explain this dif- 
ference. However, as shown in the bottom plot of Figure 9, the 
B-I color distributions of the two sets do agree. Furthermore, 
the R-I color distributions do not show any inconsistency be- 
tween the two sets. K-S tests indicate that the chance that the 
I-H-0,5 color distributions of the cluster galaxies and the caU- 
brators are consistent is less than half of those for other colors. 

Virgo and Ursa Major Clusters: Another test that indi- 
cates the preference for the 1-band results is to examine RIH_o.5 
TF distances to the Ursa-Major and Virgo clusters. Pierce and 
TuUy (1988) published RI magnitudes of nearly 20 galaxies in 
each of Ursa-Major and Virgo clusters. Also H-0.5 magnitudes 
are available from the AHM survey. We applied our calibrated 
BRI TF relations to these samples to determine the distances. 
The i/-o.5 magnitudes for these galaxies were taken from the 
AHM survey. For both clusters, we find that the //_o.5 distances 
are systematically larger, by 0.15 mag, or 7-8% in linear dis- 
tance. 

From these comparisons, we tentatively conclude that the H- 
band TF distances tend to be overestimated and should be used 
with caution: that is, there is something systematically different 
about the H-0.5 magnitudes of the nearby calibrators and distant 
cluster galaxies. The I-band distances are consistent with those 
of B-, V- and R-band results, and the H-band Hq is inconsis- 
tent with other Hq values by 2-4 a. Although they may be less 
affected by extinction, the //_o.5 magnitudes are dependent on 
isophotal diameter measurements and the concomitant aspect 
ratio corrections. There may be a systematic difference in //-0.5 
magnitudes of local caUbrators and cluster galaxies. This will 
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require further investigation, but is beyond the scope of this pa- 
per, and the discrepancy is Hkely not solved without a consistent 
set of total H-band magnitudes. However, the measurement of 
total H magnitudes is made difficult by a sky background which 
is more than 100 times brighter at H than I, and needs to be cor- 
rectly subtracted from the low surface brightness profiles in the 
outer parts of galaxies. Until better photometry than the current 
H database becomes available, the photometric systematics we 
have struck here may set the limits on Hq measurements from 
the TuUy-Fisher study. 

5. DISCUSSION 

We have derived the value of Hq by focusing on the survey 
by G97 who presented the 21cm linewidth and 1-band photom- 
etry data for several hundred galaxies in clusters out to ^10,000 
km s"' . In the previous section, this value was then compared 
with other surveys of different wavelengths. In this Section, we 
discuss how biases and systematic uncertainties can affect the 
value of Ho, concentrating on the 1-band database. 

5.1. Biases 

5.1.1. Consistency between the Calibration and Distant 
Galaxies Samples 

Not only do the local calibrators and cluster sample galaxies 
have to be on the same photometric and kinematical systems, 
they also need to be of similar type. The B- / color as a func- 
tion of linewidth has often been used as this is a good indicator 
of the variations in the recent star formation rate (~1 Gyr). It 
is also sensitive to the internal extinction correction. Unfor- 
tunately for majority of the galaxies in the G97 I-band survey, 
there are no B magnitudes available. We thus compare our local 
calibrators with galaxies in the Virgo and Ursa Major clusters 
whose photometric and Unewidth data have been published by 
Pierce & TuUy (1988). These clusters should be representative 
of other clusters of galaxies that are used in the two all-sky sur- 
veys we are calibrating, especially since the spiral galaxies used 
are not located in the immediate centers of the clusters, but in 
the "infall" regions. 

In Figure 10 (bottom), we showed the B—I color distribu- 
tions of Virgo and Ursa-Major cluster galaxies as a function 
of linewidth. The local calibrators were over-plotted by solid 
circles. The two cluster samples - Virgo and Ursa Major - 
show no significant difference between their galaxy popula- 
tions. Furthermore, we see no significant difference between 
the local calibrator sample and the two clusters. This was 
not the case when Pierce & TuUy (1988) first examined this 
diagram; they reported that six local calibrator galaxies from 
their sample were found systematically on the blue edge of this 
color-Unewidth relation. The enlarged sample here indicates 
that there is no systematic difference between the local calibra- 
tors and cluster galaxy samples. 

Another important issue that needs to be addressed is the 
difference in environment between the calibration and cluster 
galaxies. Of 21 calibrators, only seven are members of well de- 
fined groups/clusters of galaxies: the Virgo, Fornax, and Leo I 
clusters. The rest are members of loose groups or are field 
galaxies. The application of the TF calibration based on such a 
sample to the distant cluster data could potentially introduce an- 
other systematic uncertainty, as cluster and field galaxies could 
be inherently different. Using these seven calibrators only and 
fixing the slope to that obtained from the cluster data, we de- 
rived an I-band TF zero point (for a 50% linewidth relation) of 



-21.44±0.17 mag. This agrees within Icrofthe zero point ob- 
tained in Equation 11, suggesting that within our limited sam- 
ple, there is no significant evidence for environmental depen- 
dence in TF relations. A more rigorous, statistically better de- 
fined analysis is presented by Giovanelli (1997b) who exam- 
ined the 1-band TF residuals for galaxies grouped by the pro- 
jected radial distance from the center of clusters. They found 
no zero point offset between subsamples. In light of these re- 
sults, we make no attempt to separate the field galaxies from 
the cluster galaxies in the cafibration. 

5.1.2. Selection Biases in the I-band 

In measuring the Hubble constant using distant cluster data, 
it is important to account correctly for selection effects which 
include the Malmquist bias and cluster population incomplete- 
ness bias. The former arises in field surveys, for example, 
where galaxies are distributed uniformly and thus there are 
more galaxies per unit volume at larger distances. If each true 
distance is scattered randomly, then more galaxies are scattered 
"downward" (i.e. smaller distances) than "upward". Thus, the 
distances are underestimated. This type of bias was first exam- 
ined in detail by Lynden-Bell et al. (1988). Later, several au- 
thors expanded on this subject, including Gould (1993), Landy 
& Szalay (1992) and Willick (1994). As for the cluster pop- 
ulation incompleteness bias (CPIB), if not properly corrected 
for, the clusters would appear brighter than they actually are 
because only the brightest galaxies in them are observed in a 
flux-limited sample, leading to a biased value of the Hubble 
constant. 

There are numerous theoretical and empirical studies that 
have dealt with the CPIB. Teerikorpi, Bottinelli and collab- 
orators introduced a quantitative approach to this probelm 
(Teerikorpi 1984, 1987, Bottinelli et al. 1988). Willick (1994) 
also presented a more extensive quantitative analysis of the sta- 
tistical bias. We will show in this section using the formalism 
presented by G97 that the CPIB in our sample is small, affect- 
ing the derived value of the Hq by less than 5%. 

Federspiel et al. (1994, hereafter FST94) investigated the 
bias properties of the TF distances using a field sample of 1355 
galaxies from the Mathewson-Ford-Buchhom survey. One of 
the main conclusions presented by that study was that in a flux- 
limited sample, at a given redshift (which is to first order, dis- 
tance), the Hubble constant is multi-valued for different val- 
ues of rotational velocity and apparent magnitude. The trend is 
that for galaxies with similar rotational velocities, the Hubble 
constant rises as the redshift becomes larger, and for a given 
redshift, a higher Hq is found for smaller galaxies (smaller ro- 
tational velocities) (for details, see Figures 1 and 2 of FST94). 
They also showed that the zero point of the TF relation is sen- 
sitive to the redshift, and to apparent magnitude as well. Such 
biases clearly need to be corrected for if field samples are to be 
used to measure the Hubble constant. 

We show, however, that for our cluster samples of the I-band 
spiral galaxies, these biases are small. Here, we examine the 
two distant cluster surveys which we are calibrating to measure 
Ho, to see to what extent they may be affected by the incom- 
pleteness biases discussed by FST94. 

In Figure 11, I-band TF relations are shown for the G97 
cluster sample. Absolute magnitudes were calculated assum- 
ing //o = 100 km s"' Mpc"' and that velocities (in the CMB 
reference frame) give distances to first order. The choice of Ho 
here is not critical, as we are interested in examining the rela- 
tive offsets of different subsamples. In the left panel of Figure 
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10, subgroups are divided by their velocities, and the mean of 
each bin is shown. With the exception of the nearest sample 
(cz < 2000 km s"'), there is no significant difference among the 
subsamples. The slight offset in the zero point of the nearest 
sample is most likely due to its complex velocity field not be- 
ing characterized effectively by the CMB reference frame. The 
right panel of Figure 10 shows subgroups divided by apparent 
magnitudes. If the I-band TP sample is being severely affected 
by its flux limit, then we should observe a systematic trend in 
the zero points of the subsamples. That is not the case, however 
The four subsamples, although some suffer from small number 
statistics, show no significant dependence of TF zero point on 
the apparent magnitude. 

Lack of selection biases suggests that the incompleteness 
bias in the G97 1-band sample is not severe, especially if the 
TF sample is restricted, for example to those galaxies with 
larger rotational velocities (log W > 2.35 - see Figure 7). In the 
next subsection (5.1.3), we will show quantitatively that clus- 
ter population incompleteness bias for our cluster samples is at 
most 0.03-0.05 mag, or 1-2% in distance, and hence in the 
value of Hq. 

5.1.3. Cluster Population Incompleteness Bias 

We showed in the previous section that the selection bias with 
respect to apparent magnitude or redshift is negligible. In this 
section, we take a more quantitative approach to correct for se- 
lection biases. 

There are several approaches to correct for CPIB. Giovanelli 
et al. (1997b) presented a thorough, comprehensive method to 
correct for this bias, which we apply to the I-band survey. The 
basic idea is that the observed luminosity function for each clus- 
ter is compared with the expected one derived from a theoretical 
luminosity function such as the Schechter function. Then a ran- 
dom sample is drawn from the model to represent the observed 
cluster. After 1000 iterations, the zero point of the TF relation 
obtained from the simulated sample is compared to that of the 
observed one, to estimate the bias. The detailed description of 
this CPIB method is found in Giovanelli et al. (1997b). 

The cluster population incompleteness bias amounts to, on 
the average, ~ 0.03 - 0.05 mag, corresponding to ~ 1 - 2% 
in distance. Unlike what some authors have claimed in the 
past (e.g. Sandage, Tammann and collaborators), incomplete- 
ness bias is almost negligible. For example, Sandage (1994) 
had suggested that one of the reasons that the TF method gave 
a high Hubble constant value was because only the brightest 
galaxies in the cluster were observed: in particular, he believed 
that only the "upper" half of the TF relation was being consid- 
ered. This is not what is actually happening. These simulations 
show that the galaxies with log W > 2.5 are virtually unaffected 
by the cluster population incompleteness bias. Thus the overall 
bias per cluster amounts to only a couple of percent in distance. 

5.2. Systematic Uncertainties in //(> 

It is important to address how systematic uncertainties can 
affect Hq. In this section, we focus on the following: (1) 
the slope and shape of the TF relation, concentrating on how 
much a use of a wrong slope value can affect the Hq, and also 
whether the application of quadratic fits would significantly 
change the Hq values; (2) the metallicity dependence of the 
Cepheid period-luminosity relation; (3) the internal extinction 
correction method used; and (4) velocity flow field models, fo- 
cusing on how dependent the value of Hq is on the velocity 



reference frame adopted. Table 5 summarizes all the values of 
Ho derived for each test, details of which are described below. 

5.2.1. Dependence ofHo on TF Slope and Shape 

The effect of slope uncertainties on the value of Ho is small. 
We estimated the Hubble constant with three different assump- 
tions for the 1-band slope of -9.00,-10.00, and -11.00. We 
obtain 75 ± 4, 73 ± 2, and 71 ± 2 km Mpc"^ respectively. 
Errors quoted here include only the random uncertainties. In 
this exercise, we changed the value of the slope by 1.0; in re- 
aUty, this corresponds to more than IOct, which suggests that a 
use of incorrect TF slope will affect Hq insignificantly, perhaps 
2-3% at most. 

In Appendix C, quadratic TF relations were derived for each 
wavelength. To demonstrate that the changes in the shape of 
the TF relation affects Hq very Uttle, we applied the quadratic 
equation to the I-band cluster data. Note that the equation (C4) 
was derived for 20% linewidth. Here, we used a similar relation 
derived for the 50% linewidth. This yielded Hq = 73 ± 2 km 
s"' Mpc"^ 

5.2.2. Effects of Metallicity Dependence of the Cepheid PL 

Relation on Ha 

Whether the Cepheid period-luminosity relation depends on 
the metallicity of the Cepheid variable stars has been a partic- 
ularly important concern for the Hq HST Key Project. It may 
be one of the larger sources of systematic errors in the Cepheid 
distance scale. The metalhcities of the Key Project galaxies 
vary almost by an order of magnitude. All the distances have 
been determined with respect to that of the LMC, for which we 
adopt (m-M)() = 18.50±0.13 mag. However, the metal abun- 
dances of the TF calibrators cover a wide range, up to five times 
higher than LMC. Several observational and theoretical studies 
claim to have constrained this metalhcity effect (see Kennicutt 
et al. 1998). Some estimates suggested a dependence of as 
much as 0.8 magnitude per factor of 10 increase in abundance 
(e.g. Gould 1994). Because most of the galaxies used in our TF 
calibration are more metal rich than the LMC, if such a steep 
dependence indeed existed, it could potentially modify the TF 
zero points significantly, and hence the derived value of Hq. 

We have re-derived TF relations using the metallicity de- 
pendence of 7 = 6im-M)o/6[0/H] = -0.24 ± 0.16 mag/dex, 
derived from the differential test of Cepheid distances to two 
fields in M 101 (Kennicutt et al. 1998). Figure 12 shows the 
1-band TF relation without corrections for the metallicity de- 
pendence of Cepheid PL relation (solid circles) and with correc- 
tions (open circles). The metallicity-corrected set of BVRIH-0 5 
TF relations are: 

55,2 = -8.02(±0.72)(logWfo-2-5)-19.78(±0.11) (12) 

= -9.00(±0.83)(logW^o-2-5)-20.34(±0.12) (13) 

2 = -8 .89(±0.72)(logW5) - 2.5) - 20.67(±0. 11) (14) 

/^2 = -9.42(±0.76)(logW2o-2.5)-21.15(±0.12) (15) 

i/fo.5,z = -ll-21(±0.87)(logW2'i,-2.5)-21.80(±0.14) (16) 

We have also re-derived the 1-band TF relation using the sec- 
ond calibration method which was actually applied to the dis- 
tant clusters to derive Hq. This method incorporated the cluster 
galaxy data, using the 50% linewidths: 

7^2 = -10.00(±0.12)(logW'5'o-2.5)-21.41 (17) 
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Table 5 

Systematic Errors affecting tlie I-band Hq 



Errors 



Errors 



Ho 



TF Slope 

-11.00 71 ±2 

-10.00 73 ±2 

-9.00 75 ± 2 

Quadratic Fits 73 ± 2 

With Metallicity Correction 70 ± 2 
Without Metallicity Correction 73 ± 2 



Exlinclion Correction 

T98 73 ±2 

GiovaneUi et al. 74 ± 2 

Han 73 ± 1 



Velocity Models 
Local Group 
CMB 

Infall Model 



71 ±3 
73 ±2 
74±3 



'Unless otherwise indicated, the velocities in the CMB reference frame were 
used. 



The metallicity-corrected zero points above are brighter, as 
most of the galaxies in this particular sample are more metal 
rich than the LMC. Following the same exercise shown in Sec- 
tion 3.2, we have calculated the TF zero point by adopting the 
slope measured from the cluster galaxies. The value of Hq de- 
rived using the Cepheid distances corrected for their metallicity 
dependence is Hq = 70 ± 2 km s"' Mpc~', ~ 4% smaller than 
one obtained without any corrections. 

5.2.3. Effects of the Uncertainties in the Adopted Internal 

Extinction Corrections on Hq 

One of the largest corrections applied to the raw magnitudes 
is the internal extinction correction, for which we adopted the 
expressions derived by T98. The internal extinction for highly 
inclined galaxies amounts to more than 0.5 mag. In Table Al 
(in Appendix A), we listed values of internal extinction for each 
Cepheid caUbrator galaxy determined via three separate meth- 
ods: Tully et al. (1998), GiovaneUi et al. (1994) and Han 
(1992). The large dispersion between these methods suggests 
that the internal extinction corrections are not being determined 
uniquely. Can the choice of extinction correction method affect 
the measured value of HqI 

To quantify such an effect, we derived Hq using all three 
above methods of internal extinction correction. For each set of 
calculations, magnitudes of both caUbrator and cluster galaxies 
are treated using the same corrections for consistency, then the 
slope and zero point of the 1-band TF relation are estimated, as 
is the incompleteness bias for each cluster, and finally Hq is cal- 
culated from the measurements of TF distances to the clusters. 
For three extinction correction methods, T98, GiovaneUi's and 
Han's, we obtain Hq = 73 ± 2, 74 ± 2 and 73 ± 1 km s"i Mpc"i 
respectively. As we had expected, Hq is virtually insensitive to 
the extinction correction method used. 

Furthermore, there is of course the inclination uncertainty 
which propagates through the extinction correction eventually 
to the error in Hq. A full error analysis is provided in the next 
section. 

5.2.4. Dependence o/Hq on the Large-Scale Velocity Field 

One of the major remaining uncertainties in the determina- 
tion of the Ho is the correction of the observed velocities of 
galaxy clusters for large-scale peculiar motions. It is now clear 
that there are motions on scales of tens of Mpc's with ampU- 
tudes up to few hundreds of km s~^ with respect to the CMB. 



However, the exact nature of these motions is stiU unclear 
(Lauer and Postman 1994; Riess, Press and Kirshner 1995), 
especially in terms of what precisely are the causes of these 
motions. In Table 3, we listed three velocities: with respect to 
the Local Group reference frame, CMB frame, and the mod- 
eled flow field, all of which were used to derive Hq for each 
cluster for I-band surveys.The observed heUocentric velocity 
(Vh) of the galaxies are corrected to the centroid of the Local 
Group following Yahil, Tammann and Sandage (1977), and also 
to the CMB frame via the relation from GiovaneUi et al. (1998). 
The flow field model adopted in this paper is a multi-attractor 
model which included the Virgo cluster, the Great Attractor and 
the Shapley concentration. The details of the flow field model 
are presented in Mould et al. (1999). 

The flow field model is explained briefly as follows. It is a 
multi-attractor linear model based on Han and Mould (1990) 
and Han (1992). The model is comprised of three attractors: 
(1) the Virgo cluster at (RA,DEC,Vh) = (12.472/z, +12°.67, 
1035km s~'), with an assumed radius of 10 degrees; (2) the 
Great Attractor at (13.333/;, -44°00, 4600km s"'), with an as- 
sumed radius of 10 degrees; and (3) the Shapley concentration 
at (13.50/1, -31°0, 13800km s"'), with a radius of 12 degrees. 
The flows toward each attractor are assumed to be independent 
of each other, thus the total velocity correction is a vector sum 
of all three pecuUar motions. 

As shown in Table 5 in which we have listed Hq estimates 
using three different velocities, Hq is not sensitive to the veloc- 
ity field adopted. The value of Ho remains unchanged whether 
we adopt the velocities in the CMB frame of reference or those 
derived using the infaU model. 



6. UNCERTAINTY IN THE HUBBLE CONSTANT 

Understanding all the sources of errors, both systematic and 
random, and how they propagate into the uncertainty in the 
value of the Hubble constant is almost as challenging as the 
actual measurement of Hq itself. Table 6 shows how the un- 
certainty in Ho is derived. Here, we focus on the derivation of 
uncertainties for the 1-band survey, as that for the other surveys 
should be similar. 

The first part describes the errors in the Cepheid distance 
scale. Much of this has been discussed in detail in previous 
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papers of this series. The largest systematic error originates in 
the uncertainty of the distance to LMC, with respect to which 
all of the Cepheid distances have been measured. Combining 
the LMC distance uncertainty with other errors such as those 
in the HST photometric zero points, extinction treatment and in 
the photometry of Cepheid variables themselves, the total ran- 
dom and systematic errors of the Cepheid distance scale add up 
to 0.08 and 0.16 mag respectively. 

Next, the uncertainties in the TF distance scale are consid- 
ered. Using simulations, we have determined that the errors in 
the zero points and slopes of the TF relations are 0.09 and 0.13 
mag respectively. For the H-band TF relation, however, we as- 
sign a larger error for its zero point (0.25 mag) because of the 
systematic uncertainties in the photometry of galaxies between 
the calibrator and cluster samples. This was discussed in Sec- 
tion 4.2. A typical error in the measured apparent magnitude of 
cluster galaxies is -^0.04 mag. We have also assigned a typical 
error of 15 km s~' to the linewidth measurements, which is a 
conservative estimate. For a typical galaxy in our sample with 
logW ~ 2.6, a 15 km s"' error is equivalent to 0.04 mag. There 
is in addition an uncertainty of 0.04 mag in the internal ex- 
tinction correction, which was estimated from comparing three 
different correction methods (Appendix A). Furthermore, the 
Galactic reddening adds an uncertainty of 0.02 mag. These un- 
certainties, however, are reduced by the number of target galax- 
ies in each cluster, y^. Thus the total random error in the TF 
magnitudes per cluster is only 0.03 mag, assuming a typical 
number of galaxies per cluster of 16. The total linewidth error 
is 0.01 mag, again assuming n = 16. 

We learned in Section 3.4 that the intrinsic dispersion of the 
TF relation is ~ 0.20-0.25 mag. Again, the uncertainty in the 
TF cluster distance due to the intrinsic dispersion is reduced by 
a factor ^Jn. The total random error in the TF distance moduli 
is therefore the photometric error, the linewidth error multipUed 
by the slope, and the intrinsic dispersion combined in quadra- 
ture. Note that we have excluded the error in the TF slope. We 
have shown in Section 5.2.1 that the slope uncertainty has a 
negligible effect on the value of Hq. 

The total systematic error in the TF distance modulus con- 
sists of three terms. The first is the error in the TF zero point 
calibration, which we estimated to be 0.09 mag using Monte- 
Carlo simulations. Another term is the systematic uncertainty 
arising from the internal extinction correction. We have ex- 
cluded any galaxies that required an internal extinction correc- 
tion of 0.6 mag or more, however, there still might be a slight 
systematic error due to the fact that the mean extinction cor- 
rections are not the same for the clusters and calibrators. We 
have included this uncertainty (0.04 mag) in the total system- 
atic error of the TF distance moduli. Adding these two terms 
in quadrature with the systematic error of the Cepheid distance 
scale (0.16 mag), we estimate the total systematic error of the 
TF distance scale to be 0.22 mag. 

Finally, the uncertainty in the value of Hq can be determined. 
In addition to the total random error on the TF distance moduli 
derived above, the random error of the cluster velocities needs 
to be considered. In the last section, it was shown that the 
value of Ho was affected by the velocity flow field adopted, due 
to uncertainties arising from not understanding fully what the 
sources of peculiar motions of clusters are, and also due to the 
non-uniform distribution of clusters with respect to the CMB 
dipole. Giovanelli et al. (1998) report that the distribution of 
pecuUar velocities of clusters in the CMB reference frame has 
a dispersion of 300 km s"^ . We assign this 300 km s~' as an 



additional random uncertainty in the determination of the Ho 
value. 

Combining all the errors, we obtain Ho = 73 ± 2 ± 7 in the 
I-band (the first and second uncertainties correspond to ran- 
dom and systematic respectively), 72 ± 2 ± 1 1 in the B-band, 
68 ± 2 ± 10 in the V-band, and 67 ± 2 ± 10 in the H-band. 
Again, larger systematic errors for the H-band result are due 
to its less certain TF relation zero point. We have also adopted 
larger systematic errors for the B- and V-band results as for 
these, the photometric data for calibrator and cluster galaxies 
were not derived in a consistent manner, although both use to- 
tal magnitudes. We take the weighted average of these four 
estimates to derive the value of the Hubble constant, where the 
weights are defined as the inverse square of the quadrature sum 
of the errors that are intrinsic to the TF relations. These include 
the random error on the TF moduli (R2.4 in Table 6), the sys- 
tematic error on the TF zero point (82. 1), and the systematic 
error in the internal extinction correction (S2.2). In Table 6, the 
error estimates for the I-band TF relation were listed. For other 
wavelengths, the extinction correction uncertainties are scaled 
by 2.75 for B-band, 1.5 for V-band and 0.5 for H-band. Thus, 
the total random error on the TF moduli are 0.12 mag, 0.11 
mag, 0.12 mag and 0.12 mag for B, V, I and H-band TF rela- 
tions respectively. As indicated above, the systematic error on 
the TF zero point is 0.13 mag for the I-band, while for others, 
we have adopted 0.25 mag. The weighted value of the Hubble 
constant is therefore: 

7/o = 71 ±4 (random) ±7 (systematic) km s"' Mpc"'. (18) 

7. SUMMARY 

The HST Key Project on Extragalactic Distance Scale has 
contributed enormously in building the cosmic distance scale. 
Taking the Tully-Fisher relation alone as an example, the num- 
ber of its calibrators have quadrupled as a result of Cepheid 
observations using the HST. This significantly improved the ac- 
curacy of the TF relation zero point, and effectively the value 
of Hq. Furthermore, numerous all-sky surveys, such as the one 
we used in this paper by GiovanelU et al., allowed us to ex- 
tend the galaxy distances far out enough such that their peculiar 
velocities are significantiy small compared to their recession 
velocities. 

In this paper, we have derived the value of the Hubble con- 
stant using the I-band TuUy-Fisher relation. The BVRIH_o.5 
TF relations were calibrated using ground-based photometry 
data from Macri et al. (2000) and Aaronson, Huchra and Mould 
(1982,1986), and Cepheid distances published as part of the se- 
ries in Papers I through XXIII of this series. These were then 
applied to the published cluster surveys of Giovanelli et al. (I- 
band), Bothun et al. (B-band and V-band), and Aaronson et al. 
(H-band) to derive the value of Ho. The clusters extend out to 
Vcmb — 10,000 km s~'. Taking the weighted average of Hub- 
ble constants derived using four different surveys, we obtained 
Hq = 1\ ±4 (random) ±7 (systematic) km s"' Mpc"'. 

Table 7 lists published Ho values derived using the TF rela- 
tion. Focusing on the values published during past five years, 
the estimates vary from the low 50s to the low 80s. There are 
mainly three reasons for this variation. First is the calibration 
of the TF relation. Mould et al. (1996) showed that doubling 
the number of calibrators changed the zero point of the TF rela- 
tion by nearly 0.3 mag, consequently lowering the value of Ho 
by 15%. Another example arises when comparing our I-band 
Hq value with that of Giovanelli (1997). As shown in Table 5, 
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Table 6 
Error Budget (I-band) 





Source 


Errof^ 


Notes 


1. 


Gcphcid Distance Senile 








A. LAIC True Modulus 


±0.13 


Adopted from Westerlund (1997) 




B. LMC PL Zero Point 


±0.02 


From \f adore & Freedman (1991) 


Sl.l 


T Mr* PT Sv^tpinnrir TTrrnr 


±0.13 


A nnH R jiHHf»H in nnnHmtiirp 




D. HST V-Bcmd Zero Poinf 


±0.03*' 






E. HST I-Band Zero Point'' 


±0.03*' 


^JEi^{\ -Rf +E^R^, R = A(V)/E(V-I) = 2.45 


SI. 2 


Systematic Error in the Photometry 


±0.09 


Rl.l 


Random Error in the Photometry 


±0.05 


From DoPHOT/ALLFRAME comparison 




E Rv Differences Between Galaxy and LMC 


±0.014 


See Ferrarese et al. 1998 for details 




G. Errors In the adopted value for Ry 


±0.01 


See Ferrarese et al. 1998 for details 


R1.2 


Random Error in the Extinction Treatment 


±0.02 


F and G added in quadrature 




H. PL Eit (V) 


±0.05' 




/. PL Fit (I) 


±0.04' 




R1.3 


Random Error in the Cepheid True Modtxlus'' 


±0.06' 


H and I partially correlated 




Total Random Error 


±0.08 


Rl.l, R1.2 and R1.3 added in quadrature 


SpL 


Total Systematic Error 


±0.16 


Sl.l and S 1 .2 added in quadrature 


2. 


TuUy-Fisher Distance Scale 








A. Error in the measured apparent magnitude 


±o.io7V«' 






B. Error in the corrected linewidth 


±0.047 






C. Error in the internal extinction correction 


±0.047 


average value from Table Al 




D. Error in the Galactic reddening 


±0.02 


Typical value given in §6 


R2.1 


Random Error in the TF Magnitudes 


±0.03 


A,C and D added in quadrature with « = 16 


R2.2 


Random Error in the linewidth 


±0.01 


n assiuned to be 16 


R2.3 


Intrinsic Dispersion of TF Relations 


±0.20/ 


^jRl.l^ + iW.O X «2.2)2 +fi2.3^ 


R2.4 


Random Error on the TF moduli 


±0.12 


S2.1 


Systematic Error on the TF Zero Point 


±0.13 


from Equation 9 


S2.2 


Systematic Error in the internal extinction correction 


±0.04 


from Table Al multipUed by AA/ (cluster - calibrators) 


S2.3 


Systematic Error due to the metallicity dependence of 








Cepheid distance scale 


±0.08 


see Section 5.2.2 


Rtf 


Total Random Error on the TF Distance Moduli 


±0.12 


R2.4 


Stf 


Total Systematic Error on the TF Distance Moduli 


±0.22 


S2.1, S2.2, S2.3, and Spl added in quadratore 


3. 


Hubble Constant 






R3.1 


Random Error on the TF Cluster Velocities 


±300'' 


in km/sec 


Rho 


Total Random Error on Hq 


±2 


^y{R3.l Ivf + (OAmoRTF flVN, in km/sec« 


Sh(| 


Total Systematic Error on Hq 


±7 


O.A(>H(iStf in km/sec/Mpc 



"Errors are in magnitudes for sections 1 and 2, and as indicated for section 3. 

''Contributing uncertainties bom the Holtzman et al. (1995) zero points, and the long-versus short uncertainty, combined in quadrature. 

'The values quoted are typical but individual cases vary slightly. 

"'The partially correlated nature of the derived PL width uncertainties is taken into account by the (correlated) de-reddening procedi 
coupled with the largely 'degenrate-with-reddening' positioning of individual Cepheids within the instability strip. 

'GiovaneUi's analysis of the velocity field. 

'Each cluster is assumed to have n TF galaxies. 

= 5000km/s, and N clusters of galaxies are used to derive the Hq. 
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Table 7 



Estimates of Hq using TuIIy-Fisher Relation 


Ho 


Method" 


Sample 


Source 


50±4 


BTF 


Virgo 


Sandage & Tammann 1976 


80±8 


BTF 


Virgo 


Tully & Fisher, 1977 


76±8 


BTF 


Virgo 


BottineUi & Gouguenheim 1977 


61 ±4 


IRTF 


Virgo and Ursa Major 


Aaronson. Mould & Huchra, 1979 


65 ±4 


IRTF 


Virgo 


Mould. Aaronson & Huchra, 1980 


95±4 


IRTF 


4 clusters <6000 km/s 


Aaronson et al. 1980 


63 ±10 


(H-r)vs log W relation 


Virgo 


TuUy, Mould & Aaronson 1982 


71±2 


VTF.RTF.IVTF 


Clusters <6,000 km/s 


Visvanatiian, 1982 


82± 10 


BTFJRTF 


Groups <6000 km/s 


Aaronson & Mould. 1983 


74± 11 


VTF.RTF.IVTF 


Clusters <6000 km/s 


Visvanathan. 1983 


55 ±9 


IRTF 


Coma tied to Virgo 


Sandage & Tammann. 1984 


91 ±3 


IRTF 


20 Sc rs< 13,500 km/s 


Bothun et al. 1984 


92±1 


IRTF 


Clusters < 11,000 km/s 


Aaronson et al. 1986 


85± 10 


B.R.IRTF 


Virgo and Ursa Major 


Pierce & Tully 1988 


57 ± 1 


BTF, IRTF 


Virgo 


Kraan-Korteweg. Cameron & Tammaim. 1988 


68 ±8 


BTF 


Virgo 


Fouqueetal. 1990 


92!?^ 


BTF 


Coma 


Fukugita et al. 1991 


73±4 


IRTF 


Clusters < 11,000 km/s 


Mould et al. 1996 


69±8 


IRTF 


Virgo, Ursa-Major 


Shanks 1997 


82 ±10 


BTF 


Virgo 


Yasuda, Fukugita & Okamura, 1997 


69±5 


ITF 


Clusters <~10,000 km/s 


Giovanelli et al. 1997a 


53 ±5 


BTF magnitudes 


klun'' 


Theureau et al. 1997 


57 ±5 


BTF diameters 


klun'' 


Theureau et al. 1997 


57 ±7 


BTF 


Clusters < 1 1 .000 km/s tied to Virgo 


Federspiel. Tammaim & Sandage 1998 


56±3 


BTF 


KLUN (HST Calibration) 


Theureau 1998 


51±4 


BTF 


KLUN (HIPPARCOS calibration) 


Theureau 1998 


53 ±3 


BTF 


Field galaxies <5000 km/s 


Tammann 1999 


56±3 


BTF 


Clusters < 1 1 .000 km/s tied to Virgo 


Tammann 1999 


77 ±4 


ITF 


Clusters <8.000 km/s 


Tully 1999 


76±3 


ITF 


Clusters <8000 km/s 


Madore et al. 1999 


52±5 


BTF inverse diameter 


KLUN 


Ekhohnetal. 1999 


53±6 


BTF inverse magnitude 


KLUN 


Ekholm et al. 1999 



'ITF: I-band TF; IRTF: IR TF; BTF: B-band TF; VTF: V-band TF; RTF: R-band TF; IVTF: 10500A 

''KLUN stands for Kinematics of the Local Universe, which consists of a sample of 5171 spiral galaxies out to ~5000 km/s. 
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using Giovanelli et al's extinction correction method, and ap- 
plying that to all the clusters (no velocity cutoff), we obtained 
Ho = 76± 3 km s~' Mpc~', whereas Giovanelli (1997) reported 
Ho = 70 ± 5 km s"' Mpc"^ . Even though these two values agree 
with each other within Icr, the difference is attributed to a dif- 
ferent calibration. Giovanelli (1997) used a smaller sample of 
1 1 galaxies, instead of 21 that we used in this paper Using the 
calibrators of GiovanelU (1997), we do in fact retrieve the same 
Ho value. Increasing the number of calibrators again by 70% 
(from 11 to 19) contributes considerably to the change in the 
value of Hq. Furthermore, the same argument applies for the 
difference between the values of Hq obtained by Aaronson et al 
(1980) and this paper which used the same set of clusters. The 
calibration adopted by the former consisted only of two galax- 
ies, M31 and M33. The decrease in Hq is almost entirely due to 
the fact that these two galaxies he well below the TF relation, 
thereby systematically yielding a lower Ho value. 

The second reason behind the various inconsistent estimates 
is the treatment of the cluster population incompleteness bias. 
Sandage, Tammann, Federspiel and collaborators apply a large 
incompleteness bias to their B-band TF appUcation, largely 
based on the assumption that the TF intrinsic dispersion in B- 
band is as large as 0.7 mag. We have shown that the intrinsic 
TF dispersion is only ~0.30 mag, even in B-band. In fact, there 
is no anti-correlation between the wavelength and the disper- 
sion, as some papers had suggested earUer (e.g. Pierce & TuUy 
1988). 

Third, some of the B-band TF distances to clusters have been 
measured with respect to the Virgo cluster whose distance mod- 
ulus was assumed to be 3 1.60 ±0.08 mag (Tammann 1999). 
When there are nearly two dozen TF calibrators, there is no 
need to estimate TF distances based on one cluster distance. 

Results from the KLUN (Kinematics of the Local Universe) 
survey have given Ho values around 50-55 km s"' Mpc"' (cf. 
Theureau et al. 1997). This is a B-band, 21cm survey of ~2000 
spiral galaxies out to 2000 — 3000 km s~^ The latest results 
from this survey (Ekholm et al. 1999) reported Ho = 52 ± 5 km 

Mpc"' using the inverse diameter B-band TF relation, and 
53 ± 6 using the inverse magnitude B-band TF relation. The 
KLUN database includes field galaxies, which need to be cor- 
rected for Malmquist bias. They also found that their calibrator 
sample and the field sample followed different diameter slopes, 
and subsequently corrected for this effect. Both of these cor- 
rections decrease the value of Ho. Our TF analyses suggested 
that the cluster incompleteness bias was minimal, requiring a 
correction of only 2-3% in distance. Furthermore, we did not 
find a significant difference between the slopes of the calibrator 
and cluster samples. 

We have also explored how systematic uncertainties can af- 
fect estimates of Ho. The main results are sunmiarized as fol- 
lows: 

(1) Intrinsic dispersions of the TF relations are a ~ 0.19- 
0.25 mag for BVRIH-0,5. The dispersion does not seem to de- 
crease significantly with wavelength as reported earlier. 



(2) No color dependence of the TF relation was observed. 

(3) The cluster population incompleteness bias was found to 
be small, typically ~ 0.05 mag at most per cluster in distance 
modulus. 

(4) The dependence of Hq on the TF slope is negligible. 

(5) We found that our estimate of Ho is not sensitive to the 
choice of internal extinction correction. The estimates of Ho us- 
ing three different extinction correction methods all agree well 
within 1<T uncertainties. 

(6) A very small dependence of the value of Hq on the metal- 
hcity dependence of the Cepheid PL relation was observed. 
When a metallicity dependence of 7 = 0.24 ±0.16 mag dex~^ is 
adopted, Hq decreases by 3%. 

(7) Ho is not sensitive to the velocity field model adopted. We 
have used the CMB reference frame as a default case in mea- 
suring Hq. However, we have also examined the case in which 
the flow field was estimated using a three-attractor model. That 
yielded an Hq that is exactly the same as that obtained using the 
CMB velocities. 

(8) The larger two databases, based on the I-band and H- 
band surveys, which are hkely more statistically reUable than 
the other two, give Hubble constants that are different by 10%. 
We suggest that this difference arises from systematic photo- 
metric differences between the calibrators and cluster galaxies 
in H_o.5 magnitudes, which are circular aperture magnitudes. 
This clearly needs to be revisited in the future, possibly by ob- 
taining deep H-band total magnitudes consistently for both the 
calibrators and cluster galaxies. 
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APPENDIX 

A. INTERNAL EXTINCTION CORRECTIONS 

We briefly review three different methods for correcting for internal extinction. The "internal extinction" correction here strictly 
refers to the relative correction of spiral galaxy magnitudes to the face-on orientation; we do not apply an absolute, total extinction 
correction to any of our galaxies. 

The first such correction to consider is that of Han (1992). Using a sample of 284 spiral galaxies in I-band, he presented a set of 
parametric expressions as follows: 

r -0.731og(fe/fl) T<3 
Afnri=\ -0.901og(fe/fl) T = 4,5 (Al) 
[ -0.5Uog(b/a) T>6 

where (b/d) is the minor-to-major axis ratio, expressed by (1 - e). 

G97 examined again several hundred Sbc-Sc galaxies and derived a more complicated set of parametric formulae for internal 
extinction corrections which depended on both the inclination and the size of the galaxies. The basic idea behind their correction is 
that the light traverses a longer path in a larger galaxy, increasing the internal extinction; the larger the galaxy is, the more correction 
is needed. The internal extinction expression given by G97 is: 

A;^^] = 7(logH^)log(fl//7) (A2) 

where 7 is parametrized as a function of linewidth, whose values are given in Table 3 of G97 and vary between 0.57 and 1 .70. In 
addition to the size-dependent correction, G97 suggests a constant offset for galaxies of earlier types. 

C 0.10 T=l,2 

j3 = I 0.32 T = 3 (A3) 
[ 0.0 otherwise. 

The total extinction is thus given by + (}. 

B. INCLINATIONS 

In Table Bl, we summarize the inclination information compiled from previously published data for each galaxy. The columns 
include: (1) name of the galaxy; (2) photometric inclination angles published with references indicated by numbers in brackets; and 
(3) kinematical inclination angles. 

C. QUADRATIC TULLY-FISHER RELATIONS 

Several papers have investigated the possibility of a Tully-Fisher relation being a quadratic relation (cf: Mould, Han and Bothun 
1989). Pierce and TuUy (1988) noted that the TF relation at faint end tends to curve downward. Giovanelli et al. (1997) however 
noted that their sample of several hundred galaxies showed no such trend. We have fitted a quadratic equation to each BVRIH 
Tully-Fisher correlation and obtained the following results where w = logWjoyc 

21 -19.8380 -7.94259 8.11011 17 -20.4028 -8.41765 6.44754 21 -20.6891 -9.34836 9.29448 21 -21.1989 -10.19832 11.5500 20 
-21.8610-12.1775 12.7621 



=-19.84-7.94w + 8.11w^ 


(CI) 


=-20.40 -8.42w + 6.45w^ 


(C2) 


7?^ =-20.69-9.35w + 9.29w^ 


(C3) 


7^ = -21 .20 - 10.20W -t- 1 1 .55w2 


(C4) 


/f5o.5 = -21.86- 12. 18w+ 12.76w2 


(C5) 



In Figure 13, these fits are shown by solid lines. In comparison, linear fits (Equations 6-10) are drawn by dashed lines. 

If one were to use the quadratic equation instead, distances to individual galaxies may change slightly, especially for those with 
large linewidth (log W > 2.6) or smaller linewidths (log W < 2.3). However, its effect on the value of Hq should be negligible for two 
reasons: (1) the distance to the cluster should not change in the mean, and (2) we are restricting the cluster sample galaxies to those 
in the range of 2.3 < logl^ < 2.6. Because the number statistics are not sufficient to conclude whether the TF relation is quadratic 
or not, and also because the calibrator sample spans a rather small range of Unewidths, we will not use above equations (CI) through 
(C5) in our further analysis. 
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Table B 1 
Inclination Angles 



NGC 


Photometric'* 


Kinematical^ 




(Degrees) 


(Degrees) 


NGC 224 


72 (i,5), 78 (2,3,4) 


75-78 (18) 


NGC 598 


51(7), 55(5), 54(2) 


55(6) 


NGC 925 


55 (7), 52 (8,9), 51 (70), 54 (77), 57 (2) 


55 (72) 


NGC 1365 


57 (1.13), 61 (14), 63 (77), 44 (75) 


40(76) 


NGC 1425 


62 (7), 66 (75), 64 (77) 




NGC 2090 


59 (7), 60 (13), 62 (77) 




NGC 2403 


53 (7), 55 (3), 60 (4), 61 (2) 


60 (72), 50 (6) 


NGC 2541 


58(7) 




NGC 3031 


59 (1,9) 


58(6) 


NGC 3198 


65 (1,5,3) 


70 (72) 


NGC 3319 


56 (7,79) 


59(26) 


NGC 3351 


47 (7,20) 


34 (6), 40 (27) 


NGC 3368 


45(7) 




NGC 3621 


50 (1,13), 62 (3), 65 (77) 


65 (22) 


NGC 3627 


62(7) 


65(6) 


NGC 4414 


54 (7) 


50-70 (23) 


NGC 4535 


43 (7), 45 (5) 


40 (24) 


NGC 4536 


64 (5), 65 (7), 69 (20) 




NGC 4548 


36 (7) 


42 (24) 


NGC 4725 


45(7) 


53 (72) 


NGC 7331 


70 (7.5), 71 (7,23) 


75(6) 



" (1) BottineUi et al. 1984, (2) Shanks et al. 1997, (3) Huchtmeier & Richter 
1988, (4) Rood & WiUiams 1993, (5) Biviano et al. 1990, (6) Garcia-Gomez 
& Athanoussla 1991, (7) Krumin & Salpeter 1979, (8) Davis et al. 1983, (9) 
Giraud 1986, (10) Giovanelli & Haynes 1985, (11) TuUy 1988, (12) Wevers 
et al. 1984, (13) Reif et al. 1982, (14) Schoniger & Sofue, 1994, (15) Bureau 
et al. 1996, (16) Jorsater & van Moorsel 1995, (17) Mathewson et al. 1996, 
(18) Corradi & CapaccioU 1991, (19) Davis & Seaquist 1983, (20) Jackson et 
al. 1989, (21) Buta 1988, (22) Walsh 1999, (23) Braine & Combes 1993, (24) 
Guhathakurta et al. 1988, (25) Dumke et al. 1995. (26) Moore & Gottesman, 
1998 



Sakai et al. 



21 




2.4 2.6 2.8 2.4 2.6 2.8 2.4 2.6 2.8 

log W (20%) log W (20%) log W (20%) 




2.4 2.6 2.8 2.4 2.6 2.8 

log W (20%) log W (20%) 



Fig. 1 . — BVRIH_o 5 Tully-Fisher Relations for spiral galaxies with Cepheid distances, using 20% linewidth. Solid lines represent the bivariate fits, while the 
dotted and dashed lines represent inverse and direct tits. 
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Fig. 2. — Observed minus predicted absolute magnitude as a function of color. Although the sample is small, there is no obvious dependence of the TF relation 

on color. 
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Fig. 3. — Color dependence on inclination for TF calibrator galaxies. Solid circles represent the correlations after correcting for the internal extinction using 
TuUy et al.'s method (top), Han's method (middle) and GiovaneUi et al.'s method (bottom). Open circles are for the colors uncorrected for internal extinction. If the 
internal extinction correction were "correct", this correlation should have zero slope. 
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log W (20%) log W (20%) 



Fig. 4. — (left): B, I and H-band Tully-Fisher relations in which the photometric inclination angles were used. The barred galaxies are indicated by open circles. 
The solid lines represent the TF fit through all the galaxies, while the dashed lines the fit using only the barred galaxies, (right): B, I and H-band TF relations in 
which the kinematical inclination angles were adopted for barred galaxies only. Their TF fits are shown by solid lines. The dotted lines are the fits to all the galaxies 
using the photometric inclination angles (solid lines on the left plots). The zero point change in, for example, the I-band TF relation is 0.01 mag. 
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Fig. 5. — The V-H and B-H color distributions of the local calibrators. The barred galaxies are indicated by open circles. On the left, the magnitudes were 
corrected using the photometric inclination angles, while for those shown on the right side, the kinematical inclination angles were substituted for the barred 
galaxies. 
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Fig. 6. — I-band Tully-Fisher relation for the cluster survey galaxies. The solid circles represent those galaxies classified as barred. 
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Fig. 7. — Examples of comparisons of the calibration galaxies (solid circles) with the cluster galaxies in the I-band survey (crosses). 
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Fig. 8. — Comparison of distributions of internal extinction corrections: using Giovanelli et al (top), Tully et al. (middle) and Han's method (bottom). The soUd 
histograms represent the distributions for cluster galaxies, while the dotted-dashed histograms are for the caUbrators. In our analysis, we exclude any galaxies whose 
internal extinction exceeds 0.6 mag. 
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Fig. 9. — A log Z)- log V relation for I-band clusters (top). A straight line represents Ho = 73 km s"' Mpc"' . The bottom figure shows Ho as a function of velocity 
in the CMB reference frame. Taking the mean of clusters with Vcmi, > 2000 km s"' , we obtain Hq = 73 ± 2 (random) km s" Mpc~' . 
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Fig. 10. — Top: Comparison of distance moduli derived from the 1-band and H-band TF relations for galaxies that are included in both surveys. The solid line 
represents an one-to-one correspondence and is not a fit through the data. Middle: The I-H^j color distribution as a function of logarithmic linewidth for Virgo 
and Ursa Major galaxies (open squares) from Pierce & Tully (1988), and A 1367, A400 and Coma galaxies whose photometric data (open triangles) were compiled 
from the G97 and B85 databases. The local calibrators are overplotted with solid circles. Note that the local calibrators are located on the red edge of the correlation 
for the cluster galaxies. Bottom: B—I color as a function of linewidth for the same set of galaxies as in the middle figure. The local calibrators are overplotted (solid 
circles). Four galaxies used by PT88 whose Cepheid distances were measured from ground are indicated by solid circles overplotted with open circles. 
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Fig. 1 1 . — I-band Tully-Fisher relations for samples grouped by velocities in the CMB reference frame (left), and grouped by apparent magnitudes (right). The 
absolute magnitudes were calculated assuming Hq = 10()km sec"' Mpc"' and that velocities give distances to first order. With an exception of the nearest sample, 
there seems to be no dependence of Tully-Fisher zero point on the sample. For the nearest sample (< 2000kms"' ), the offset is likely due to the fact that its velocity 
field is not characterized well by the CMB model. 




Fig. 12. — I-band-log W20 Tully-Fisher relation (solid circles). Open circles present the luminosity-linewidth relation for same galaxies, except that their 
distances have been adjusted to incorporate the metallicity dependence of the Cepheid PL relation, using 7 = 0.24 mag dex"' (see text). 
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Fig. 13. — The solid line in each Tully-Fi.sher relation repre.sent.s a quadratic fit while the dashed line is the linear relation (Equations 6 - 10). Applying the 
I-band quadratic TF relation to the cluster data, we obtain Hg = 75 ± 5 km s"' Mpc^' . 



